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Abstract 



"Naturalness of electroweak physics within minimal supergravity" 

A thesis submitted for the degree of Doctor of Philosophy 

Sebastian Cassel, St. John's College 

Trinity 2010 



Low energy supersymmetry is motivated by its use as a solution to the 
hierarchy problem of the electroweak scale. Having motivated this model 
with naturalness arguments, it is then necessary to check whether the ex- 
perimentally allowed parameter space permits realisations of the model 
with low fine tuning. The scope of this thesis is a study of naturalness 
of the electroweak physics in the minimal supergravity model. The latest 
experimental constraints are applied, and the fine tuning is quantitatively 
evaluated for a scan across the parameter space. The fine tuning of the 
electroweak scale is evaluated at 2-loop order, and the fine tuning of the 
neutralino dark matter thermal relic energy density is also determined. 
The natural regions of the parameter space are identified and the associ- 
ated phenomenology relevant for detection discussed. Naturalness limits 
are also found for the parameter space and spectrum. The minimum fine 
tuning found is 1 part in 9 when dark matter constraints are neglected, 
and 1 part in 15 when dark matter constraints are satisfied. For both 
cases, the minimum fine tuning is found for a Higgs mass of 115 GeV 
irrespective of whether the Higgs mass constraint is applied or not. The 
most natural spectrum includes light superpartner fermions, and heavy 
superpartner scalars. Minimal supergravity currently remains viable with 
respect to naturalness and a natural realisation may be discovered within 
the next couple of years. 
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Chapter 1 



Introduction 



A theory of nature must account for the experimentally observed symmetries. These 
include the relativistic invariance of spacetime along with an independent set of sym- 
metries, SU{3) X SU{2) X [/(I), of the fundamental matter. The Standard Model 
(SM) of particle physics has proved extremely successful in quantitatively matching 
experimental measurements, based on these local symmetries alone and the minimal 
particle content required to match observations. 

There are however many problems with the Standard Model that motivate de- 
velopment of 'beyond the Standard Model' (BSM) physics. An obvious flaw is the 
absence of gravitational interactions. Standard Model physics is based on quantum 
field theory, for which a consistent physical description of the non-renormalisable in- 
teractions and dynamics of gravity requires an infinite number of parameters. The 
unpredictability at ultraviolet energy scales is then interpreted as an indication that 
the model is incomplete or defective. The Standard Model also fails to explain neu- 
trino masses or the existence of dark matter. 

Hints for the existence of further BSM physics may be derived from the many 
apparently coincidental features of the Standard Model. For example, the SM gauge 
couplings seem to all converge at a high energy scale on solving their renormalisation 



group equations. The quantum field theory anomahes of the quark sector cancel 
those of the independent lepton sector. The electric charges of all known matter is 
an integer multiple of a common quantised unit, leading to the electric charge of 
the electron to be exactly opposite that of the proton. Questions such as why there 
aren't more (or less) fundamental particles in nature, or why there are precisely three 
families of matter with identical sets of quantum numbers are also not answered by 
the Standard Model. 

On examining the Standard Model, there are also various aesthetically unpleasant 
theoretical features. There are a large number of parameters which must be fixed from 
experiment. The values of these parameters vary over several orders of magnitude. 
Some parameters such as the degree of violation of the discrete charge and parity 
symmetries of the QCD sector are restricted to zero within experimental error without 
any theoretical justification. This leads to asking why the SM parameters take the 
values they do, and if there is a deeper reason to expect the required pattern of values 
to describe our universe. 

The concept of naturalness has been a driving principle in recent decades for many 
proposals of BSM physics. The physical principles governing a model such as special 
relativity in an arbitrary quantum field theory permit a vast family of theories which 
seem equally viable before experimental validation. The philosophy of natural model 
building is that further physical principles must be identified in order to restrict 
the viable theories of nature. In a "theory of everything," a desired quality is that 
the physics of our universe is fully explained with little need for calibration of its 
parameters, if any. 

A measure of naturalness can be associated with the sensitivity of observable 
quantities to changes permitted in the underlying fundamental parameters. "Un- 
natural" theories are then deemed to require large amounts of fine tuning of its 
parameters to match observations. The Standard Model suffers from such fine tuning 



problems. An alternative view to naturalness is to follow anthropic arguments and 
accept that arbitrary patterns just occur, or that a multiverse exists and conscious 
observers are necessarily only present in universes with appropriate physics leading 
to a selection bias. It is the implications of a principle of naturalness that is explored 
in this thesis. 

With the Large Hadron Collider (LHC) beginning to open a new energy frontier 
and Fermilab still pushing up the luminosity frontier, we may soon be able to identify 
the physics responsible for electroweak symmetry breaking (EWSB) and/or discover 
BSM physics. In the Standard Model, the Higgs mechanism is used to trigger EWSB. 
The existence of the Higgs boson necessary for the Higgs mechanism remains the last 
component of the Standard Model unverified by experiment. 

It emerges that from a theoretical perspective, the SM Higgs sector is responsible 
for its worst fine tuning problem, excepting the scale of the cosmological constant. 
Thus, if a naturalness principle is to be believed, identifying natural BSM theories 
with respect to the electroweak physics and its experimental signatures may lead to 
imminent discoveries of new physics beyond the Standard Model. 

The subject of this thesis is an analysis of the degree of fine tuning present in a 
class of BSM theories with minimal supersymmetry broken through a restricted set 
of gravitational interactions. The most natural realisations of the BSM model are 
identified with indications for promising experimental search strategies discussed. In 
this introductory chapter, the formalism of supersymmetry is presented and notation 
defined. Supersymmetric Lagrangians are introduced with spontaneous supersym- 
metry breaking mechanisms and the structure of minimal supergravity (mSUGRA) 
reviewed. The hierarchy problem is then introduced and a quantitative measure of 
naturalness defined. 

Chapter 2 presents the results of evaluation of the fine tuning of the electroweak 
scale for a scan of the allowed parameter space. A combination of analytic and 



numerical results are discussed, and the controlling features of fine tuning identi- 
fied. Chapter 3 presents the results of fine tuning of the thermal relic density of the 
mSUGRA dark matter candidate with the relevant dark matter constraints applied. 
These results are compared to the fine tuning of the electroweak scale with the dark 
matter constraints then added to the results of Chapter 2. Chapter 4 discusses how 
extensions beyond the mSUGRA model may generate improved naturalness, and a 
model independent analysis is presented using an effective theory framework. Finally, 
the conclusions of this study are offered in Chapter 5. 

The novel features of this work include analysis of the fine tuning of the electroweak 
scale at two loop order, a quantitative study of the effect of higher dimension operators 
on fine tuning, inclusion of the latest experimental constraints, and dependence of the 
fine tuning of electroweak physics on the composition of dark matter. The results of 
previous related work is also reviewed in the following chapters. 

1.1 Towards supersymmetry (SUSY) 

When building models of physics 'beyond the Standard Model' (BSM), it is vital to 
understand the full class of theories into which the set of SM symmetries may be 
embedded. Before introducing supersymmetry, a brief review of the Standard Model 
symmetries is given where notation and conventions are defined. 

1.1.1 Poincare symmetry 

In flat spacetime, the symmetry group associated with relativistic invariance is the 
Poincare group. For a set of points in spacetime, {x'^}, the Poincare symmetry is 
defined by the quadratic form of any spacetime interval being invariant under linear 



co-ordinate transformations. 

where the metric of flat spacetime, r] = diag(l, — 1, — 1, — 1). The general form of 
the co-ordinate transformation which satisfies this constraint is a combination of a 
Lorentz transform and translation. 



X 



'^ = A^x'^ + a^ (1.2) 



for some constant 4-vector, a^. A general Lorentz transform can be constructed 
from an 5*0^(1,3) sub-group along with time and space reflections. The 5*0^(1,3) 
group belongs to the set of real general linear 4 dimensional matrices subject to 
the constraints; A'^r^A = t], detA = +1, A°q > -|-1, where the last two constraints 
label the transformation as 'proper' and 'orthochronous' respectively. The subset 
of connected continuous Lorentz transformations are referred to as belonging to the 
restricted Lorentz group. 

The generators of the Lorentz symmetry, M , are defined by A = exp (— |cj^°"M ) 

The generators can be decomposed into boosts, Ki = — Mqj, and rotations, Jj = 
hijk^jk where i, j, k G {1, 2, 3}. 

An equivalent method of expressing the above quadratic form exists in a complex 
2-dimensional space. Introducing the following mapping operators, where {a^} are 
the Pauli matrices, 

a^ = (l2,{aj) a^ = {l„{-a,}) (1.4) 



spacetime 4-vectors can be mapped between the respective spaces according to: 



"'^ - ^TriXa'^] (1.5) 



The quadratic form in eq fll.ip is then given by det (X — Y). This remains invariant 
under the following transformation. 

X' = HXH^ for some H E SL{2, C) (1.6) 

SL{2, C) is the group of 2-dimensional complex matrices with unit determinant. 
The above mappings lead to the relation, A^^^{H) = |Tr [a'^HauH^. The groups 
S0_^{1,3) and SL{2,C) are locally homomorphic, possessing homomorphic Lie alge- 
bras and so both describe the symmetries of proper orthochronous Lorentz transfor- 
mations. 

To classify the representations of the restricted Lorentz group, a common approach 
is to complexify the Lie algebra. There exists a 1:1 correspondence between the 
representations of real and complexified Lie algebras [Ij making this technique useful. 
For the complex generators, Sf" = [J^ ± iKj] /2, the 5*^ form independent SU{2) sub- 
algebras. The complexified restricted Lie algebra is then homomorphic to SU{2) x 
SU{2) allowing its group representations to be classified by the eigenvalues, [S^ , S^^ . 

The (53^, 5*3") representation has a total spin of S^ + S^ and dimension {2S^ + 
1).(2S'3' -|- 1). Parity interchanges the eigenvalues, 5*^ — ?■ S^ and is obtained via 
Hermitian conjugation. By convention, the representation (|, O) is referred to as left- 
handed, and (0, 1) as right-handed. Further examples of other representations are 
given in the following table. 

The possible dimensions of S'0_,_(1, 3) representations restrict the spin of the states 
they can describe to integer values. However, the representations of SL{2, C) also in- 
clude half-integer spin states. SL{2, C) is the spin group associated with the metric rj. 
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(0,0) 


spin-0 scalar 


V2' 2) 


spin-1 4- vector 


(1,0) 


spin-1 anti-symmetric self-dual tensor 


(0,1) 


spin-1 anti-symmetric anti-self-dual tensor 


(1,1) 


spin-2 symmetric traceless tensor 



Table 1.1: Representations of the restricted Lorentz group 

The representations of a spin group are known as spinors. The fundamental spinors 
of SL{2, C) are Weyl-type and can be used to construct higher spin states. 



(1,0)0(0,1) = (1,1) 
(1,0)® (1,0) = (0,0)0(1,0) 



:i.7) 
:i.8) 



The (1, 1) state is equivalent to the fundamental representation of S0^{1, 3). To 
construct a parity invariant state with the same spin when S^ ^ S^ , a direct sum of 
opposite parity states is required. For example, a Dirac spinor has the representation 
(1,0)0(0,1). 

The invariance under translations is included to extend the Lorentz symmetry to a 
Poincare symmetry. This leads to a second Casimir operator of the symmetry group. 
In addition to spin, the rest mass associated with a Poincare representation is then a 
necessary label. 



1.1.2 Fundamental spinors of 5'L(2, C) 

There exists four different transformation rules for fundamental representations of 
the SL{2, C) group. For H G SL{2, C), these are given by; 



self representation -0^ = Hj^ ip 

dual representation '0'"^ = (^H~^^) il)^ 



B 
B 



• complex conjugate self representation ^p'■^ = (H*)^ ip- 

• complex conjugate dual representation ijj'^ = (if*^^^) ^ip^ 

where the different representations are distinguished by upper and lower indices of 
dotted and undotted type. The indices can take values of {1, 2}. 

However, the self and dual representations of SL(2, C) are equivalent in the sense 
that there exists a matrix, e G SL{2, C) such that the transformation matrices are 
related by, eHe^^ = H"^^ . Up to a sign convention, this matrix is given by. 



^AB = (e^"") ' = -^^2 

e^^ = {e^sf = ^^2 (1-9) 

where Ug is the second Pauli matrix. The e matrix acts like a metric to raise or lower 
indices of the SL{2, C) space, ip"^ = e^^^ip^, with the adopted notation. 

Similary the complex conjugate self and dual representations are equivalent, where 



T^AB _ I ^AB^ 



e = I e I = ia2 

-^AB = {^''''y = -^^2 (1-10) 



giving, ^^ = e^^ip^. 

The complex conjugation operation does not generate equivalent representations. 
By considering a general transformation on a Dirac spinor (composed of both the 
direct sum of a left and a right-handed Weyl spinor) and its Hermitian conjugate, it 
follows that the relation, ip^ = ijj*^ {a^) must hold [2] . The representation tp is 
referred to as a left-handed state, (|,0), and tjj a right-handed state, (O, |). 



Objects without free SL(2, C) indices are Lorentz invariant. Where indices are 
suppressed, a summation convention is imphed such that ip4> = ip^4>A, and ipcp = 
il^jyct)^- These definitions remove a sign ambiguity. It has been previously determined 
that these spinor fields have half-integer spin. By the spin-statistics theorem, it is then 
necessary to consider the spinor fields as Grassman variables leading to the relation, 
ilj^cpj^ = —(t)j(ip^ {= (p'^ipj^ , and similarly for contraction of fields with dotted indices. 

The Lorentz generators of SL{2, C) in the spinor representations are for left- and 
right-handed states respectively: 






(J^^- = i{a^a'' -a^'o^) (1.11) 

To identify the type of SL{2, C) indices required for the a^ and a^ matrices, 
eqs fll.Sp and (11.61) can be inspected, noting that the index structure of SL{2, C) 
transformation matrices obeys Hj^ and //^ = {H*) ^. It then follows that (c'^)^^ 
and ((T^) are appropriate, with one dotted and one undotted index. 

1.1.3 Haag-Kopuszanski-Sohnius theorem 

The physics of particle interactions is governed by the S-matrix, as defined by: 

out(PiP2 •••|kik2 ...)i„ = (piP2 ...|5|kik2 ...) (1.12) 

where the in and out states are asymptotically defined in time by the tensor product 
of single particle states with definite momenta. 

Coleman and Mandula [3] first proved that the most general connected symmetry 
of the S-matrix consistent with (i) relativistic invariance, (ii) the S-matrix being 
non-trivial and analytic, (iii) a unique vacuum, (iv) the presence of a finite number 



of massive particle states and (v) the Lie bracket of the symmetry algebra being a 
commutator, the group generators must decompose into a direct sum of the Poincare 
generators, P'^, M'^" and an arbitrary independent set of Lorentz scalar generators. 
The latter are known as "internal" symmetry generators as they do not mix with the 
geometric symmetry generators of spacetime. 

The Haag-Kopuszahski-Sohnius theorem |3] generalises the result of Coleman and 
Mandula by neglecting (v) above, finding the generalised symmetry with a Z2 graded 
Lie superalgebra where the odd generators that mix with the Poincare generators 
transform as fundamental spinors under the restricted Lorentz group. 

A Z2 graded algebra indicates that the generators can be decomposed into two 
groups, L = Lq (B L^ and that for generators, T- G L^, the Lie bracket satisfies 
[T^,Tj] G -^^odL+j2i- '^^^^ property leads to the naming convention of generators in 
Lq being even and those in L^ being odd. 

A Z2 graded superalgebra has a Lie bracket that satisfies the following symmetry 

[T„T,.] = i-iy-' [t^,t;\ (L13) 

and the Jacobi identity becomes: 

The Lie bracket is an anti-commutator when {i,j) = (1, 1), and a commutator other- 
wise. 

The most general symmetry, referred to as supersymmetry, of the S-matrix within 
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the given assumptions has the following algebra: 



QI.Q'b) 


— a 7"/^ 


[Q2,P,] 


= 


QlQl) 


= hA^'T' 


[Qa^p,] 


= 


Ql.Ql) 


= 25"^<^^. 


[Q2,M^''] 


= H^nAQs 


m,B,] 


= ^sfQl 


[QlBi] 


= -^{S*r'Ql 



:i.i5) 



where the central charges, Z"^, commute with all the generators and the Poincare 
algebra is a subalgebra. The upper case Latin indices refer to SL{2, C) indices, the 
H, V Greek indices to Lorentz indices, and the a, /3 run from 1 to A^ where iV identifies 
the "number" of supersymmetries present. The lower case Latin indices associated 
with the Lorentz scalar generators, 5;, run to some arbitrary number also, satisfying 
the algebra \B^,B^ = if-^Bi, for some structure constants, f^^}. 

The scalar generators associated with non-trivial Si correspond to "R-symmetries" . 
For A^ = 1 SUSY, the maximal R-symmetry is U{1) and there is a null central charge. 
This R-symmetry leads to invariance under chiral rotations of the spinor generators. 
For internal symmetries, the associated Si parameters must be zero to maintain in- 
dependence from the generalised geometric symmetry. 

The Coleman-Mandula theorem states that the even generators include a Lorentz 
vector, P'^ [in the (|, |) representation], an antisymmetric tensor, Af"^ [(1, 0) © (0, 1)] 
and a set of Lorentz scalars, Bi [(0, 0)]. From group theory arguments, it then becomes 
clear that the only odd generators that could mix with the Poincare generators belong 
to the (|,0) and (O, |) representations. 

To identify the superalgebra in eq fll.lSp as the most general form, consider the 
Lie bracket of two odd generators, Q, for an example in the required reasoning. 
The Lie bracket can only produce even generators due to the algebra grading and 
consideration of the direct product, (i, O) ® (i, O) = (0, 0) © (1, 0). By the Coleman- 
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Mandula theorem, there are no even generators in the (1, 0) representation, but there 
may exist Lorentz scalar generators restricting the central charge, Z"^ = af Bi, for 
some co-efficient, af. The Jacobi identities further restrict the central charges to be 
anti-symmetric in their indices, with the af co-efficients related to the Si parameters, 
which are separately fixed to be Hermitian. 



1.2 Representations of the SUSY algebra 

The Casimir operators of the superalgebra are given by 

P^P'' and B^B^P^P'^-B^B^P^P^ (1.16) 

where B is a generalisation of the angular momentum operator, and P is the usual 
momentum operator. The ffist Casimir operator is common to the Poincare alge- 
bra defining the rest mass of the representation, m. In the rest frame the second 
Casimir operator reduces to rn? X]fc=i ^k^k^ ^^'^ the Bf, commute with the spinorial 
generators. 

B, = \^,u,.P''MP''-\Q<y^Q (1.17) 

Bf^^f--^) = mJ,-\Q-a,Q (M 0) (1.18) 

The commutation relations with the standard angular momentum operator obey: 

[Q%-h] = U^3)aQb [QlJs] = -H^3)aQb (1-19) 

The action of the generator Qi on a Poincare state increases the spin eigenvalue, j^, 
by a half, and Q^ lowers the eigenvalue by a half (and vice versa for Q2, Q2)- Since a 
super- Poincare representation contains multiple Poincare representations of differing 
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spin, SUSY representations are commonly referred to as a supermultiplets. 

For a finite dimensional super-Poincare representation, a Poincare state within 
the supermuliplet can always be found for a given 63 value, such that when acted on 
by any left-handed spinorial operator, Q^, it is annihilated. This state is known as 
the Clifford vacuum, 1^2). The additional states with the same 63 value can then be 
found through repeated applications of the Q"^ operators: 

\m,b,;j,+ '^) ex g-|fi) 
\m,b^;J3-l) oc Q^\n) 

\m,b,;j,) ex g?gf|fi) (1.20) 

and so on. This procedure must be repeated for each 63 component in order to 
construct a complete supermultiplet. Note that in the absence of spinorial charges, 
the usual spin multiplet is recovered on iterating through the 63 values. For A^ = 1 
SUSY, the anticommutation relations imply that there are no further states for a 
given 63 in a supermultiplet beyond those listed in eq ( ll.20p . 

For massive supermultiplets, all the states in eq f ll.20p are physical as they have 
non-zero norm, allowing normalisation. Choosing the rest frame to test for non-zero 
norms, the anticommutation relations of type {Q,Q} obey {Q2,Q1} = 2m5"^5y^^. 
There are then 2^^~^ distinct complex fermionic and bosonic degrees of freedom each, 
for a, /3 G {1, ■ ■ ■ , A^} with a given 63 value. 

For massless supermultiplets where there is no rest frame, it is convenient to 
choose the frame with momentum aligned along the third spatial axis. The non-zero 
anticommutation relations of type {Q,Q} are {Qi,Qi} = AE 5°"^. The states with 
Q2 acting on the Clifford vacuum thus have zero norm and are not physical states. 
This reduces both the fermionic and bosonic degrees of freedom to 2^~^ each for 
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a given 63 value. Non-zero central charges (allowed when A^ > 1) may also further 
reduce the number of non-zero norm states present for the massive and massless cases. 
From now on, discussion will be restricted to the minimal N = 1 SUSY. This is 
phenomenologically favourable as to construct matter states in the fundamental rep- 
resentation of a gauge symmetry, Lagrangians cannot be written for supermultiplets 
with spin 1 component states. This restricts the generalised angular momentum of 
the matter supermultiplet to zero, and N < 2. The N = 2 and 63 G {0} super- 
multiplet contains opposite parity fermions, but parity must be broken to realise the 
left-handed electroweak symmetry. A^ = 1 SUSY allows for chiral supermultiplets and 
so any extended supersymmetry must break to A^ = 1 SUSY above an energy scale 
where electroweak symmetry is considered to exist. The analysis within this thesis is 
concerned with the energy range where electroweak symmetry is appropriate, and so 
A^ = 1 SUSY is now assumed. 

1.2.1 Superspace 

By introducing spinorial Grassman parameters, 6 and 6, the SUSY algebra can be 
written in terms of commutation relations alone. For N=l SUSY, 



O^Qa:Qb^^ = ^{(^^'^)P, (1-21) 



The usual technique for Lie groups of exponentiation can be used to find a represen- 
tation for the symmetry transformation operator. Ignoring Lorentz transformations, 
this operator is given by: 

L{x^,e,e) = expl-ix^'P^ + iOQ + ieQ] (1.22) 

In addition to the usual spacetime co-ordinates, this formulation suggests that the new 
Grassman co-ordinates can be considered to generalise spacetime into a superspace. 
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Using the Baker-Campbell-Hausdorff formula, e^^e^ = e ^ '^'2' ' '^"\ the action 
of two successive transformations is related to a single transformation by: 

L{0,e,e)L{x'',e,e) = L{x^ + iOa^'e - iea^'e.e + e,e + e) (1.23) 

A field in the superspace, or superfield, can be expanded by a Taylor expansion 
along the Grassman co-ordinates to determine a general form of the superfield in 
terms of component fields that only depend on spacetime: 

$(x,^,^) = (t){x) + V2e'^i)A{x) + V2e^x^{x) + eeF{x) + eeG{x) 

+ Oa^'eA^i.x) + eee'^x^^x) + eee^R^{x) + eeeeD{x) (1.24) 

where under a SUSY transformation, $' = L$L^^. This expansion is finite due to the 
symmetry properties of Grassman parameters, 0j^9^ = \^ab ^^- ^ product containing 
repeated Grassman numbers is zero. If the superfield has SL{2, C) indices, then the 
component fields simply inherit these extra indices. To transform properly under 
the full super-Poincare group, the component fields must be representations of the 
Poincare group. 

The product of an arbitrary number of superfields will also have a structure in 
terms of components that can be considered as a composite superfield. This composite 
superfield will then also transform correctly under SUSY transformations, a feature 
that will be useful for building Lagrangians. 

The action of the symmetry generators on a superfield can be reproduced using dif- 
ferential operator representations. Considering an infinitessimal transformation, the 
momentum operator, P is equivalent to —id . Similarly for the spinorial generators. 
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the following identities apply: 



Qa = -^ (^A-^</''^.) (1-25) 



Q^ = -i Id^'-iiaf'erdJ (1.26) 



where (9^ = d/dO"^ and 9^ = d/dO :. 

The effect of a SUSY transformation on the superfield, $, will be to transform the 
component fields into a linear combination of the original fields and their first and 
second spacetime derivatives. It will generally only be the 99 99 term which remains 
invariant up to a total spacetime derivative term. 

It is also useful to define covariant derivatives with respect to SUSY transforma- 
tions, where by construction {Dj^^Q^} = {D^,Q^} = {D:,Qb} = {D^,Q^} = 0. 

Da = d^ + ta^^J^d^ (1.27) 

D^ = d^ + ti(jf'9)^ d^ (1.28) 

These covariant derivatives anticommute with the symmetry transformation oper- 
ators, and so superfields that transform correctly under symmetry operations can 
be constructed from the general superfield in eq 01.241) by repeated actions of these 
covariant derivatives. 

Using the notation that D^ = D^Dj^ and D^ = D:D^ and similarly for the 
Grassman parameters, it follows that D^^D^^ = D,D^^ = 0. This property can 
be used to introduce constrained superfields, ^^ = D^ $ and $^ = D^ $, where 
D $^ = and D $^ = 0. These constrained superfields are referred to as left and 
right-handed respectively by convention, or collectively as chiral. The application of 
constraints is important to identify the irreducible superfield representations as the 
degrees of freedom present in the general superfield do not match those identified in 
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A left-handed chiral superfield must generally be a function of y^ ^ x^ + iOa^O 
and 9 to satisfy the constraint, Z) $^ = 0. 

^dv.O) = My) + ^0My) + 99F^{y) (1.29) 

^Lix,9,9) = (t)L{x)-i9a^'9d^(t)L{x)-\99 99d''d^(PL{x) 

+ ^9ijL{x) + ^99 d^ijL^y''^ + d9 Fl{x) (1.30) 

For the right-handed superfields, y^ = x^ — i9a'^9 is a convenient parameterisation 
giving, ^j^{y,9) = (j)ji{y) + V^9'ijjj^{y) + 99Fj^{y). A left handed superfield can be 
transformed into a right-handed superfield by Hermitian conjugation. 

The 9999 component is always a total spacetime derivative for a chiral superfield, 
but the 99 and 99 components also now remain invariant under SUSY transformations 
up to a total spacetime derivative. Products of superfields with equivalent chirality 
obey the same chirality constraint. A product of superfields with differing chirality 
does not obey either constraint, but the composite can still be considered as a general 
superfield for symmetry purposes. 

Note that the degrees of freedom within the superfield have been reduced by the 
constraint. However, there still appears to be too many bosonic degrees of free- 



dom in comparison with the discussion in Section 11.21 It will be determined in 
the following subsection that a supersymmetric Lagrangian restricts the F field to 
be non-dynamical. This field does remain vital though to maintain supersymmet- 
ric invariance and also to account for the difference in degrees of freedom when the 
dynamical fields are off-shell. 
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1.3 Supersymmetric Lagrangians 

Supersymmetry transformations change the component fields in the superspace co- 
ordinates. For physics to be invariant with respect to supersymmetry, the action must 
be invariant. This demands that any Lagrangian in terms of the component fields of 
a SUSY representation is invariant up to a total spacetime derivative, which vanishes 
for appropriate boundary conditions. 

As observed in the previous subsection, this can be achieved by extracting the 
9999 component of a general superfield or the 99 (99) component of a left (right) 
handed chiral superfield, or any product of such fields. Extraction of the components 
is possible using integration over the Grassman variables, where by definition. 



d^9 [f{x) + 9^ gA{x) + 99 h{x)] = h{x) (1.31) 

f (f9\f{x) + 9^g^{x)+99h{x)\ = h{x) (1.32) 

d^9[f{x)+ ■■■ +9999h{x)\ = h{x) (1.33) 

for arbitrary functions f{x),g{x),h{x). 

As a first step to a physical Lagrangian, consider the product of a left-handed chiral 
superfield and its Hermitian conjugate. Using eq (11.301) and Grassman identities, we 
find 

Jd'9i^^y^^ = -ct>ld,d^cl>L-ii,^a^d^ij^ + FlF^ (1.34) 

which contains the standard kinetic terms for a complex scalar field, 0, and a Weyl 
fermion, ipj^. 

We can now associate mass dimensions with the fields and parameters. As the 
action is dimensionless in "natural" units (where h = c = 1), and the d^x measure 
has mass dimension —4 (with d^ having mass dimension 1), the scalar field, 0, and 
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chiral superfield, $, must both have mass dimension 1. The Weyl fermion, ip, has mass 
dimension | and the non-dynamical field, F, has mass dimension 2. We can also assign 
a mass dimension to the Grassman parameters with 6 and 6 having mass dimension 
— I and the Grassman measure having opposite mass dimension for eqs fll.3ip to 

firaaj) to hold. 

For a renormalisable model with a Hermitian Hamiltonian, a supersymmetric 
Lagrangian must then take the form, 

f d^xd^e\i<^Ly<^L+ (^^W^('^L) + h-C-)] (1-35) 

where W, known as the superpotential, is a polynomial in left-handed chiral super- 
fields up to 0(<l'|) as it has mass dimension 3. Note that there can be no terms 
with derivatives of the component fields of a chiral superfield generated by Grassman 
integration of the superpotential, given eq f ll.SOp . It will be shown that the superpo- 
tential is responsible for the interactions and mass of the chiral superfield component 
fields, hence the name. 

To include gauge interactions and dynamics, it is necessary to find another type 
of constrained superfield with the appropriate degrees of freedom. For this, consider 
applying a reality condition to a general superfield, V = V\ restricting the vector 
component field, A^, in eq f ll.24p to be real. V is referred to as a vector superfield 
as it will contain a physical vector component but not because the superfield has any 
Lorentz indices (which it does not). 

A superfield generalisation of a gauge transformation, or "supergauge transforma- 
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tion" , can be proposed to obey the following rules. 



e^' = e-'^\^ e'^ (1.36) 

e-^' = e-^'^e'^e*^' (1.37) 

^'l = e-^^-^L (1-38) 

^'r = ^Re'"-' (1.39) 



where A is some arbitrary left-handed chiral superfield. From the Baker- Campbell- 
Hausdorff formula, the vector superfield transforms as follows. 

V = V + i{A-A^) + ^[V,A + A^+0 (A^) (1.40) 

Note that V still satisfies the reality condition. The notation above implicitly assumes 
summation over the gauge group generators, T". In more explicit notation, let V = 
2gV"'T"' and A = 2gA°'T°'. Using the notation that the scalar component of A" is ^", 
the transformation of the A component follows: 

K = A^ + d^ (Re [2^]) - 9 r"' A\ x Re [2^] + O(e^) (1.41) 

which does indeed give the usual gauge transformation for the vector component field, 
where f°^'^ is the structure constant of the gauge algebra. The degrees of freedom in A 
can be chosen to leave Re [^] arbitrary, but set others to cancel many of the component 
fields in V . This choice of supergauge is known as the Wess-Zumino gauge. In this 
gauge the non-zero components of the vector superfield are: 

Y^^ = ea^'eAi{x) + eee\''{x) + eeex'{x) + eeeeD''{x) (1.42) 

The action of a SUSY transformation will generally change the supergauge, but this 
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supergauge choice is useful for identifying the number of degrees of freedom in the 
constrained superfield. 

To construct a superfield that will allow for kinetic terms of the component fields, 
consider the following definitions, 

W^ ^ -i^V^D^e^ (1.43) 

W^ = -{D^e^D^e-^ (1.44) 

where under a supergauge transformation, W^ = e~*^ Wj^ e^^. Note also that due 
to the action of the contracted covariant derivatives, the W superfields obey the 
constraints D W^ = and D W^ = 0. The W^ (W^) are thus left- (right-) handed 
chiral superfields with a SL(2, C) index. Their 66 and 66 terms will be invariant 
under SUSY transformation up to a total spacetime derivative. 

To identify the components of these chiral superfields, perform an expansion of 
the exponentials in eqs (I1.43P and (ll.44p . 

e-^D^e^ = D^V+'^[D^V,V]+0{V^) (1.45) 

For Abelian symmetries, the commutator is zero. In the Wess-Zumino gauge, the 
0{V^) corrections are zero, which gives the relation, 

W^{y,6) = X\{y) + 26^D\y)-{a^'^6)^F-^M 

+ i66 a^^ [d.X'^^y) - g r'^ A^X"" {v)] (1-46) 

To maintain Lorentz invariance, there must be no free SL{2, C) indices in the 
Lagrangian. It follows that the only renormalisable term in the Lagrangian involving 
these chiral superfields must be composed of W^^W^ and its Hermitian conjugate. 

For A^ to be identified as a gauge boson, the vector superfield must have mass 
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dimension 0, permitting the exponentiation of the superfield for use in Lagrangians. 
Similary, the Wj^, W^ must both have mass dimension |. 

A supersymmetric Lagrangian can now be written with dynamical matter and 
gauge fields. For an arbitrary set of matter chiral superfields, <l>j, the Lagrangian 
density is given by: 



£ = /"ci4^(^$Je2^^"^"<l>, + 2fi:\/"+ [^2j^i [l + ^] Py"^W^^ + W^($)) +h.c. 



F:F^ - ti>^ a^ D^^, - ^\ D^D^cj^, + 2g D^ (0* T>,) - V2g {^XT^^ + h.c.) 



+ 2D'^D" - iX" a^d^y - \ F;^F"'''^ + ^FI^^^F"'"'' + 2k D" 



+ 



dWi(j)) p _ 1 dW((j>) 



^i^j + h.c. 



;i.47) 



where 9 is not a Grassman variable, but a CP violating parameter for the gauge 
field dynamics. The use of Latin indices for the different chiral superfields' F-terms, 
Fj, distinguish them from the gauge vector field strength, F" , which possess Greek 
indices. Similarly, the usual spacetime covariant derivative, D is separate to the 
vector superfield's D-term, D". 

The non-dynamical fields, F^ and D", can be integrated out of the Lagrangian 
using their on-shell equations of motion. This proceedure generates the following 
Lagrangian which describes equivalent physics. 



- V2g {i'r\''T'^<P, + h.c.) - i (g (0* T'^^J + nf 



a rpa fAU 



dW((P) 



\2 d^^dcPj ^t'f^J ^ ^'-'- 



:iA8) 



In the presence of multiple gauge fields, the exponent of the exponential is replaced 
by a sum over the gauge superfields. Note also that only f/(l) gauge symmetries will 
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permit a V"" term in the Lagrangian. A general superpotential will be of the form 

W{<!>) = h% + i /i'^$,<l',. + ^ Y'^^<^^<^^<!>f^ (1.49) 

For a single matter chiral superfield the interactions present in eq (I1.48P include 

™|' = /iV + >^/i0' + x<^' + /^(/^ + 2/i0 + F0') (1.50) 

^^^.^. = H^'^^ + Y^^<l>) (1-51) 

where the scalar, 0, and Weyl fermion, ip, have degenerate masses and interactions 
controlled by only one further parameter. The mass degeneracy for degrees of freedom 
within a supermultiplet was expected from the previous arguments on SUSY algebra 
possessing the Casimir operator, P^P^. If this were not the case, the supersymmetry 
must have necessarily been broken. 

In a general A^ = 1 SUSY theory, a f/(l) R-symmetry is permitted. This involves 
chiral rotations of the Grassman parameters. The Grassman measure, d^9^ is invariant 
under such an operation, as are $^$ terms. However, for invariance under the R- 
symmetry, the superpotential must have an R-charge of magnitude 2. R-charges must 
then be assigned to the matter chiral superfields and only terms with appropriate 
composite R-charge are allowed in the superpotential. 

Noether currents can be associated with continuous symmetries of Lagrangians, 
and the derivation holds for invariance under transformations with a spinorial gen- 
erator. Consider a component field transforming as 0' = (p + e^ (^0)a ^'-'^ some 
infinitessimally small Grassman parameter, e"^, leaving the Lagrangian invariant up 
to some total spacetime derivative, C = C + ^^Q^Ja for some JT^. A conserved 
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current is formed by the definition, 



fA = J^A-yiT^rfi-^i^^dA (1-52) 



'A Z^ a(a^^^^ 



wliere 5^j^ = and the sum is over all component fields. Similarly, an independent 
Noether current can be defined for transformations, 0' = + e : (A0) . The SUSY 
Noether currents, or "sup er current s" , have three SL{2, C) indices (as a Lorentz vector 
index, /x, requires two SL{2, C) indices), and have spin | each. 

1.4 Supersymmetry breaking 

The presence of supersymmetry forces a mass degeneracy between pairs of bosonic 
and fermionic degrees of freedom. As experimental limits exclude this possibility, then 
if SUSY is a symmetry of our universe it must be spontaneously broken in the phase 
that we observe. A sign of supersymmetry breaking is that the vacuum transforms 
non-trivially under SUSY transformations. For an infinitessimal transformation, 

|0)' = |0) + ^(6^Q^ + e^g-^)|0) (1.53) 

and so for preservation of supersymmetry, Q^ and Q^ must annihilate the vacuum. 
As the Hamiltonian can be written in terms of these spinorial generators using the 
SUSY algebra, 

(0|i^|0) = (0|P°|0) = 1{0\{{Q^,Q.} + {Q^,Q.})\0) (1.54) 

the energy of a vacuum invariant under (global) SUSY transformations must be zero. 
A theory with supersymmetry does not allow arbitrary shifts in the energy scale to 
set the vacuum to have zero energy. 
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The terms in the potential of the Lagrangian in eq f ll.48p which only include 
scalar degrees of freedom are collectively known as the scalar potential. These can be 
expressed solely in terms of the auxiliary fields, which are non-dynamical. 

Vs = J]|F/ + 2D"D" (1.55) 

i 

As the gauge superfield D-terms are real, each term in the scalar potential is positive 
definite. For the vacuum to have non-zero energy and thus SUSY to be broken, one of 
the auxiliary fields must have a non-zero vacuum expectation value (vev). The terms 
in the potential with non-scalar fields must necessarily have null vevs to maintain 
Lorentz invariance. 

To break supersymmetry whilst maintaining internal symmetries (in the presence 
of non-trivial interactions), multiple matter chiral superfields are required. To realise 
this, consider the case of a single matter superfield. If this solitary matter field is 
charged under a gauge group, quantum anomalies necessarily break the internal gauge 
symmetry, irrespective of whether supersymmetry is broken. The matter field is then 
excluded from being charged under the gauge group. If the matter field has mass 
and/or self interactions, the F-term vev is then always null, as the complex scalar 
field will acquire a vev in order to minimise the scalar potential. 

The requirement for multiple matter fields is not a problem phenomenologically, 
but to illustrate SUSY breaking in a toy model, this must be realised. For an exam- 
ple of F-term SUSY breaking, consider the O'Raifearteigh model with the following 
superpotential 

W = m ^a'^g + A $1 (<l>^ - /i^) (1.56) 
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and canonical diagonal kinetic terms. The scalar potential is given by, 



2AV' 



A V' + k i<Pl <P3 



2 

+ (2Am 0^0203 + h.c.) + 4A2 |0i03|' + A^ |03 



V 



/ 



+ m^ |02| 



V</^3/ 

4 



:i.57) 



It is then impossible for any choice of (^3) to obtain a vacuum solution with zero 
energy. Supersymmetry is broken spontaneously as although the Lagrangian is still 
invariant under SUSY transformations, the vacuum is not. The bosonic degrees of 
freedom have their degeneracy split with respect to the fermionic degrees of freedom. 



For example, when m"^ > 2A^/i' 



) \S^2,3 



) = and the scalar mass spectrum includes 



ym^±2A^/i^ whilst the non-zero fermion masses present are m. 

To demonstrate D-term SUSY breaking, consider the Fayet-Iliopoulos model which 
has a U{1) gauge field with a non-zero k, term in eq fll.48p . There are also two mat- 
ter chiral superfields, $_|_ with respective U{1) charges, ±g. The superpotential can 
only contain the term, W = r7i$^$_. The Fayet-Iliopoulos model then contains the 
following equations of motion for the auxiliary fields. 



-m(f)-. 



D 



K + gq 



:i.58) 



Note that it is impossible to simultaneously set (F^) = (D) = for any non-zero 
K and so supersymmetry must be broken. If m"^ > gqK,/2, the vacuum is found at 
(0_i_) = which preserves the gauge symmetry and it is the D-term with a non-zero 
vev. The scalars again receive a mass splitting, A/m^±^^K/2. If m^ < gqK/2, both 
SUSY and the gauge symmetry are broken in this model. 

The splitting of bosonic degrees of freedom is a common feature in spontaneously 
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broken SUSY theories. For a tree level Lagrangian, the following sum rule is obeyed 

Supertrace = ^ [(-1)2-^(2J + 1) m^] = -A'^gT^{D'') (1.59) 

where the sum is over all superfields that may mix and J is the spin of the component 
field. 

After accounting for quantum corrections, the superpartner states of Standard 
Model fields should have been observed if the vev of an auxiliary field associated 
with a SM field is responsible for spontaneous SUSY breaking. This introduces the 
need of a hidden sector with no SM charges to be responsible for SUSY breaking in 
order to avoid the quantum corrected sum rule constraint. The breaking of super- 
symmetry would then be communicated through messenger states which couple to 
both the Standard Model and hidden sectors, and/or certainly through gravitational 
interactions. To avoid modelling the hidden sector physics, effective theories can be 
considered. 

1.5 Minimal supergravity 

Using the covariant derivative, Z), the action related to the Lagrangian in eq (11.481) 
can be equivalently expressed as follows, 

S = j (fx(fe(-^[^^e^].^^ + W{m) + \W''^Wl\+\i.c. (1.60) 

where the CP violating parameter, 6*, is dropped for simplicity, as is the k parameter 
that would lead to SUSY breaking. 

For a theory to include gravity, non-renormalisable interactions must be present. 
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The most general Lagrangian consistent with global supersymmetry is of the form: 

^global = /rf^xrf^^(-f /C([$teV^]^,$,) 

+ W ({$}) + \ /„, ({$, W%}) W'^W'X) + h.c. (1.61) 

where the /C function is real and includes terms with arbitrary positive powers of 
the chiral superfields and their Hermitian conjugates. The superpotential is now a 
polynomial of infinite order and the gauge kinetic function, /^^ is also a polynomial 
function of chiral superfields. To maintain the appropriate mass dimension of the 
action, non-renormlisable terms have co-efficients with negative mass dimension. 

It is now stated without proof (see [5]), that for a theory with local supersym- 
metry, the general form of the action invariant under general superspace coordinate 
transformations is of the form: 

^locai = j [d'xd'dE)[lMp[MpD'~m) e^v[-\K[[^^e''l,^,)] 

+ W ($,) + \ /,, ($,, ly^) W^^W'^ + h.c. (1.62) 

where K = — Slog [/C/(3Mp)]. The E superfield is introduced to form an integral 
measure which is invariant under the coordinate transformations. This is the gener- 
alisation of the ^—detg^^, factor in general relativity. The TZ superfield is also the 
generalisation of the Ricci scalar in general relativity. 

The Hilbert action for a spin 2 field is obtained from the zeroth order term in 
the expansion of the exponential, and so the respective component field is identified 
as the graviton. The Rarita-Schwinger action for a spin | field is also found, with 
the field referred to as the gravitino being the superpartner of the graviton. For this 
reason, the presence of a local supersymmetry leads to the identification of theories 
with such symmetry as supergravity. 
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By analogy with gauge symmetries, we must expect a dynamical spin | field when 
promoting the supersymmetry from global to local. In gauge theories, an A^j^ term 
is introduced into the Lagrangian to maintain invariance up to a total spacetime 
derivative, where j'^ is the Noether current of the symmetry. For supersymmetry, the 



Noether current was found to be spin | in Section [L3] and so a spin | field is necessary 
to couple to the current. When the gauge symmetry is spontaneously broken, the 
field coupling to the Noether current acquires a mass, and similarly the gravitino 
acquires a mass when SUSY is broken. 

In order to identify a low energy effective theory that may be derived from super- 
gravity, a minimal model is now considered. The Kahler potential, K and superpo- 
tential are chosen to be of the form below with diagonal kinetic terms as a first step 
in avoiding fiavour mixing 

i 

+ Mp' Y. ^'' (^-' ^n) ^^^j + O {M-/) (1.63) 

W = # (S) + i fi'^ (S) <l>,<l>^. + If Y'^'' (S) $,$,$fc (1.64) 

where Mp is the Planck mass, $j are some observable fields, and S„ are some hidden 
sector fields. The higher order terms in $ are ignored as they will be suppressed 
by the Planck scale and so be negligible at low energy scales. However, with SUSY 
breaking occurring in the hidden sector the same approximation cannot be made for 
a small set of non-renormalisable interactions involving S. 
The scalar potential of supergravity is given by 

V = -F, G) F' -3 Ml, e^^/^^^ + i D" (Re/ab)"' D^ (1.65) 
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where the F^ and D"' terms, and G function are given by 



G 



M, e-^/(^^i^) (G-^)J G^ + i/:,, {G-r, r A^ - (G^^); G^' ^^^, - ^ G,. ^^^ 



Qi ^rj.a- 



{^^ 



Ml 



K-\oi 



\wi{j:}{m' 



;i.66) 



using the notation that G* = 9G(0, 0*)/90j, and Gj = dG{(j),(j)*)/d(j)*. The A" are 
the gaugino fields and the component fields of $j include the scalar, 0^-, and Weyl 
fermion, ■?/'•. 

On F-term SUSY breaking triggered by one of the hidden fields, the following 
terms are generated in the Lagrangian to leading order in Mp^: 



4oft = - {m')\ <t^:<t^, - \\ M,A"A« + \ &^ 0,0^. + i, A^^' 0,0^.0, + h.c.l (1.67) 



which when canonically normalised are given by: 






\{ReQ-\F,d'Q 



5i 



m 



3/2 



+ VJMI, - {F^F"" a"a„ log K^'^) 



Y^ZAF^d"" K + \ogY^^'-\og 



J{(.i) J{U) J{(k) 



) ( 



f^i,y.,{Fmd"'{K + \ogfi'^-\og 



K^i)K^3) 



- "^3/2) ( 



(1.68) 
(1.69) 
1.70) 
1.71) 



with the caveat that O {Z^^) terms have been omitted in B^^ , and the constraint 
fab ~ fa^ab ^^ required to maintain gauge invariance of the observable fields. The 
gravitino mass, "nig ^ and Vq parameter is given by 



m 



3/2 



^ -(G>/(2M|.) 



V, = -M|e-<^>/^^(G*(G-^)jG^. + 3M|,) 



^1.72) 
:i.73) 
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The choice of eqs fll.63p and f ll.64p has thus generated a certain pattern of terms 
in the effective Lagrangian that exphcitly break supersynimetry. The K^^^ parameters 
give the curvature of the metric in field space. If a common curvature is assumed for 
all observable fields, then universal soft scalar mass terms are generated. 

The experimental constraints on flavour changing neutral currents places strong 
limits on the non- degeneracy of sfermion states with identical electroweak quantum 
numbers. The limits on third generation states are generally weaker though as, for 
example, the B^ — B^ mixing already has a large Standard Model contribution and 
significant theoretical errors. On running to the UV scale, it is important to obtain 
near universal soft masses amongst the generations for states with the same gauge 
quantum numbers. Further motivations will be provided in later sections for realising 
universality between other groups of states. 

The procedure to canonically normalise leads to the /x*-^ and V^^^ parameters 
present in the original superpotential of eq f ll.64p not corresponding to the "physical" 
parameters of the superpotential found after canonically normalising. The relation is 
not important for this discussion, but is mentioned for awareness. If the {d"^ logF*-'^) 
term in A^^^ is common to all observable fields, then the relation that A'-'^ is propor- 
tional to the physical Yukawa coupling, Y^^'^, is found. This also aids in eliminating 
the introduction of flavour changing processes. 

Due to the proportionality of certain parameters, it is convenient to introduce new 
notation where i?'-' = B^^ /x*-', and A^^'^ = A'^^ y^^'^, where the "physical" label is now 
always used. 

If the gauge kinetic function, fa = f+[ fa/^p ) ^) where the F-term of E acquires 
a vev, the gaugino soft mass terms, M^ will be proportional to /^. In the minimal 
scenario, a common factor is assumed leading to universal gaugino masses. 

Minimal supergravity (mSUGRA) is the model defined by assuming a structure in 
the non-renormalisable Lagrangian that leads to universal SUSY breaking parameters 
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at tree level. These assumptions can be motivated as the following section will discuss, 
but it is not certain that nature will be realised in such a manner. 

In addition to communicating SUSY breaking through the interactions of eq fll.63p 
and fll.64p . quantum anomalies associated with a local scaling of the superspace co- 
ordinates will generally break SUSY at loop level with a different pattern of contri- 
butions for the SUSY breaking parameters. If messenger states exist that couple to 
the hidden sector fields responsible for breaking SUSY and the observable SM fields, 
then further contributions to the SUSY breaking parameters will be present, again 
with a different pattern. It is assumed in mSUGRA that these contributions to the 
SUSY breaking parameters are negligible, and SUSY breaking is then parameterised 
at an appropriate energy scale by four universal parameters; rriQ, the soft scalar mass, 
m-^,2 5 the soft gaugino mass, Aq, the soft trilinear coupling, and Bq, the soft bilinear 
parameter. 

1.6 Minimal Supersymmetric Standard Model 

The minimal set of superfields required to contain component fields that include a 
family of SM matter states and the Higgs field are listed in the Table 11.21 There must 
then be three copies of the Q, U, D, L, E left-handed chiral superfields to reproduce 
the three families, and vector superfields must be present to contain the gauge fields. 
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Table 1.2: Representations of MSSM particle content 



Q^ contains the doublet of left-handed up and down quarks. The charge conjugate 
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of U^ contains the charge conjugate of the right-handed up- type quark (which is left- 
handed), and similarly for {Dj^y. The L^^ contains the doublet of left-handed electron 
and neutrino, and Ej^ is the right-handed electron. 

The superpartners of the SM states have a naming convention that the partners of 
quarks and leptons are squarks, q, and sleptons, /. The partners of the gauge bosons 
are gauginos, g, and the partners of the Higgs bosons are Higgsinos, h. 

In building the Minimal Supersymmetric Standard Model (MSSM), the only choice 
available to contain the SM matter is in chiral superfields. This is required since the 
fermion component field of a vector superfield must be in the adjoint representation 
of the gauge group for the Lagrangian of the vector component to respect gauge 
symmetry. None of the SM matter is in the adjoint representation of a gauge group. 
For the scalar Higgs, it is again a chiral superfield that is necessary for a scalar 
component field to be present. 

The SM only contains a single Higgs field, but in constructing the MSSM, a second 
Higgs field is required. This is a consequence of the fermion component field affecting 
the triangular gauge anomaly. Without having two Higgs superfields with opposite 
hypercharge, the SM gauge symmetries would be broken through quantum loops. The 
holomorphicity restriction of the superpotential to only include superfields with the 
same chirality for SUSY invariance also leads to the requirement for a second Higgs 
superfield in order to give mass to both the up and down quark sector through the 
Higgs mechanism. 

Having introduced the representations of the MSSM, a comment on gauge unifi- 
cation will now be offered (applying equally to the SM). One method of embedding 
the SM symmetries into a higher gauge group uses a SU{5) symmetry. A selection 
of SU{5) representations break down into SU{3) x SU{2) x f/(l) representations as 
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follows: 



5 = (3,l,|)©(l,2,-i) (1.74) 

10 = (3,2,i)©(3,l,-|)©(l,l,l) (1.75) 

24 = (8, 1,0) ©(1,3, 0)0(1, 1,0) ©(3,2,-1) ©(3, 2, 1) (1.76) 

It can be seen that the matter fields can be grouped such that {Q, U, E} G 10, and 
{D,L)} G 5. The Higgs fields are contained in 5 and 5, and the gauge fields in 24. 
The SM gauge symmetries could then be unified into SU{5) at some high energy 
scale. Other embeddings are also possible. 

The presence of gauge unification helps to justify some assumptions in minimal 
supergravity. Consider a non-renormalisable interaction, W^ S Wa where the chiral 
superfield S acquires an F-term vev. The S superfield must be in a gauge repre- 
sentation that is the symmetric product of two adjoint representations in order to 
maintain the unified gauge symmetry. For SU{5), 

(24 X 24)^ = 1 + 24 + 75 + 200 (1.77) 

If S is in the singlet representation, then the MSSM gauginos are forced to have 
universal soft masses at the scale of SUSY breaking. Soft masses for matter fields are 
similarly constrained when in a common multiplet of the unified gauge group. The 
assumptions leading to universality in mSUGRA may then be justified by further 
new physics. In general, there are 105 new real parameters required to describe the 
MSSM with explicit soft SUSY breaking terms. This is problematic for explorations 
of the parameter space and generally doesn't account for relations generated by UV 
completions with supersymmetry restored. 

The MSSM superpotential can be constructed by writing down all interactions 
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permitted by the SU{3) x SU{2) x [/(I) gauge symmetries. 

W^MssM = ^^H^H, + Y^QUH^ + Y,QDH, + Y^LEH, + W,,, (1.78) 

W^^^ e {LH^,LLE,LQD,UDD} (1.79) 

where the Yukawa couphngs, Y^^^ are 3x3 matrices in family space. 

The interactions contained in W-g^p^ are allowed by the gauge symmetries but lead 
to phenomenological problems. This is where an R-symmetry can be helpful. A 
choice for the definition of the R-charge as (— l)3('S-i)+2S'^ where 5, (L) is the baryon 
(lepton) number and S is the spin of the component field is made. The presence of 
gaugino soft masses breaks the f/(l)^ continuous symmetry allowed in A^ = 1 SUSY 
to a discrete Z2 parity. The R-parity violating (RPV) terms of the superpotential are 
then fully contained in W^p^. If this R-parity is realised, those terms are forbidden. 

This choice of R-charges assigns a common charge to all of the component fields 
present in the SM, but an opposite charge to their superpartners. The presence of R- 
parity therefore makes the lightest supersymmetric particle (LSP) stable with respect 
to decay. If the LSP has no electric charge or colour, it becomes a candidate for dark 
matter. The Standard Model itself provides no such candidate. 

However, it is not necessary to invoke R-parity to avoid the related phenomeno- 
logical problems. Another candidate for a unified gauge group is 5*0(10) that breaks 
to SU{5) X U{1) with the SU{5) breaking as mentioned before to the SM gauge 
groups. The necessary charge assignments of the U{1) forbid the terms in Wp,pv 
The 5*0(10) unified group also has the appealing feature that a whole family of mat- 
ter is contained in its fundamental representation along with a SM gauge singlet that 
may be considered as a right handed neutrino. 

The phenomological problems alluded to from VFj^pv include violation of lepton 
number [LH^,LLE,LQD) and violation of baryon number {UDD). The UDD and 
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LQD interactions lead to tree level proton decay through squark mediation, p — )• £+7r°, 
but the proton has an observed lifetime in excess of 10'^^ years. It is expected that 
gravitational interactions cause proton decay, and even Planck mass suppressed opera- 
tors can produce too rapid proton decay. The presence of renormalisable interactions 
that lead to proton decay would therefore need extremely small coefficients or the 
squarks mediating decay must be very heavy in order to satisfy experimental bounds. 
The following section will argue that both of these scenarios are "un-natural" and so 
to maintain naturalness, R-parity will be assumed. Similar restrictions are found for 
the other RPV interactions. 

1.7 Hierarchy Problem 

The hierarchy problem of the Standard Model refers to the situation that the scale 
of electroweak symmetry breaking is 0(10^ GeV) whereas gravity has an intrinsic 
scale set by the Planck Mass, Mp ~ 0(10^^ GeV), hierarchically many orders of 
magnitude larger. As the Standard Model does not account for gravity, the model 
must be considered as an effective theory with an upper limit of validity at the energy 
scale where gravitational interactions become important. 

To understand the influence of higher scale physics on low energy phenomenology, 
it is useful to consider renormalisable models with a hierarchical spectrum. The one 
loop corrections to scalar masses include terms of the form: 

5m' oc ^\j 4^^ (1.80) 

where m- is of order the mass of the heaviest virtual particle within the loop of 
the given diagram, and A- is a coupling constant. The integral is divergent and 
must be regularised. For cutoff regularisation, the loop momentum is bounded by 
a regularisation mass scale. A, and for dimensional regularisation, evaluating the 
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integral in li (= 4 — 2e) dimensions produces a finite result. A regularisation mass 
scale is also required in dimensional regularisation to correct the mass dimensions of 
the Lagrangian parameters (eg A^- — )■ yU^'^A^). The value of the integral, /, in eg 11.801 
is then given by: 



lim 



lim 

e-S>0 



167r2 V"" ""1 ""*= m'j 



:i.8ii 



inn?- 
167r2 



ih? — rn?j log 

(i + 1 - 7 + log(47r) + log ^ + 0(e))] (1.82) 



for the cutoff and dimensional regularisation schemes respectively, where 7 is the 
Euler-Mascheroni constant. 

To construct an effective theory, the degrees of freedom at the high mass scale, 
M, are removed from the Lagrangian. The contributions that this physics generates 
must be absorbed into the effective Lagrangian parameters of the reduced theory. The 
singularities are controlled through counterterms, but for both regularisation schemes 
there are 0(M^/167r^) contributions to scalar mass squared terms. The scalar masses 
then seem unprotected from acquiring masses of order close to the high energy physics 
scale. 

In the effective theory of the Standard Model, the one loop corrections to the 
Higgs mass squared, m^, is stated here with cutoff regularisation: 

^< - 6^ (3/ + g'^ + 8A - 8/^2) A2 (1.83) 

where g {g') is the electroweak (hypercharge) gauge coupling, A is the quartic Higgs 
coupling [(A/4)/i^ G £], and h^ is the Yukawa coupling of the top quark. The cutoff 
here is restricted to be less than the Planck scale, as otherwise the effective theory 
framework breaks down. 

The Higgs mass at tree level is given by, m\ = 2Af^, where the Higs vev, v = 
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246 GeV. In order to obtain a physical Higgs mass around the electroweak scale with 
new physics appearing at the Planck scale, the tree level mass squared then needs 
to be fine tuned to 1 part in ~ 10^" for cancellation of the loop corrections. This 
degree of tuning, controlled by the hierarchy of the electroweak and Planck scales, is 
considered extremely "un-natural." 

The reasoning for expecting the Higgs mass to be around the electroweak scale 
relies on demanding perturbativity and renormalisability below the Planck scale, and 
also for agreement with precision electroweak measurements. The 1-loop /3 function 
for SM Higgs quartic coupling is, 

/^r ^ sf^ (4A2 + 2/^?A - /.?) (1.84) 

For A too large, a Landau pole will be encountered before the Planck scale. For A too 
small, the top Yukawa coupling can drive A negative before the Planck scale causing 
the Higgs potential to be unbounded from below at the renormalisable level. This 
bounds the Higgs mass from above and below which for a momentum cutoff at 0{Mp) 
permits a Higgs mass between around 130 and 180 GeV [6J. As the Higgs field con- 
tributes to loop corrections of electroweak physics, there are also indirect contraints 
from precision electroweak data giving an upper mass bound of 122 (157) GeV at the 
68% (95%) confidence level [6]. 

The above arguments hold for the Standard Model being valid up to the Planck 
scale, but for BSM models where new physics is encountered below the Planck scale, 
the contraints can be relaxed to permit a SM-like Higgs mass of up to O(TeV). A 
model independent argument for new physics close to the electroweak scale can also 
be found by demanding perturbative unitarity. 

To determine the unitarity limit for a two body scattering process, it is useful to 
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consider a partial wave expansion of the cross section. 

oo 

a = if^5^(2£+l)|a/ (1.85) 



1=0 



where s is the centre of mass energy, i is the angular momentum, and the partial 
wave amplitude satisfies, \a^\ < 1. The unitarity limit is reached when \a^\ = 1. 
Using the optical theorem, a general constraint is found, |a^p = Im(a^), which given 
the structure of the partial wave amplitude implies that |Re(a^)| < 1/2. 

Considering the case of elastic scattering of longitudinal electroweak gauge bosons, 
W^, with momenta, kf^ ^ ^ (1, 0, 0, ±1), and respective polarisation vectors, e^ — )■ 
^ (±1, 0, 0, 1), it follows that the contribution from the tree level diagram with a four 
point W vertex to the i = partial wave amplitude grows as (s^/m^). Adding tree 
level diagrams with photon and Z mediation in the s- and t-channels, helps regularise 
the amplitude to grow as {s/m^r)■ For perturbativity, a sufficiently accurate result 
should be obtainable by truncating the loop expansion. However, as the centre of 
mass energy is increased, the tree level amplitude will eventually violate unitarity. 
This signals a breakdown in the viability of the perturbative analysis. 

Perturbative unitarity can be saved by a Higgs-like resonance, which modifies the 
tree level amplitude to behave as {ml/m^r). An upper limit is then placed on this 
resonance, rrif^ < O(TeV), such that the unitarity limit is not violated. The physics 
responsible for electroweak symmetry breaking is then expected from perturbativity 
arguments to be detectable below this O(TeV) scale. If perturbativity is not main- 
tained, then composite electroweak boson states will form in order to obey unitarity 
of the scattering process. This situation will not be discussed further. 

Solutions to the hierarchy problem through BSM physics are centred around mod- 
ifying the loop corrections to the Higgs mass. It is useful to note the opposing contri- 
butions from the Higgs scalar quartic coupling and top Yukawa coupling in eq f ll.83p . 
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This follows from the spin statistics introducing a relative sign difference for scalar 
and fermion loops. To avoid fine tuning of the Higgs mass, the quartic coupling 
may be very close to cancelling the contributions from the other parameters in the 
coefficient of A^. However, this scenario simply transfers the fine tuning onto the 
quartic /Yukawa coupling and will require a full loop expansion to confirm. 

Supersymmetry provides a solution which greatly improves the degree of fine tun- 
ing present. As shown in previous sections, the degrees of freedom associated with 
SUSY representations come in pairs of bosonic and fermionic states. In addition, 
the bosonic couplings are related to the fermionic couplings with the precise relation 
required to cancel the quadratic dependence on the cutoff. 

For BSM models with SUSY, the dependence of scalar masses on the cutoff scale 
is logarithmic as found for fermion masses in the SM. The SM renormalisation group 
equations are then only valid up to the energy scale where superpartner states are 
accessible. In order for the hierarchy problem to be solved, superpartner states must 
exist around the electroweak scale. The non-observation of superpartner states has 
placed lower bounds on the masses of superpartner states that are beginning to suggest 
the requirement of a small hierarchy between the electroweak scale and SUSY state 
masses. This is referred to as the little hierarchy problem. An analysis of the degree 
of fine tuning present on satisfying the experimental bounds for mSUGRA is the 
subject of this thesis. 

Alternative solutions for the hierarchy problem include positing extra compact 
spacetime dimensions. This reduces the energy scale at which gravity must become 
important, and so reduces the hierarchy between the electroweak and scale of gravity. 
Another popular model is technicolour, which suggests that electroweak symmetry 
breaking has a dynamical origin similar to the spontaneous chiral symmetry breaking 
in QCD. The electroweak scale is set by dimensional transmutation on confinement 
of the technicolour charge. Such solutions will not be discussed further. 
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1.8 Measure of naturalness 

The concept of naturalness is somewhat vague being an aesthetic quahty. A definition 
proposed by 't Hooft is [7]: 

The naturalness criterion states that one such [dimensionless] param- 
eter is allowed to be much smaller than unity only if setting it to zero 
increases the symmetry of the theory. If this does not happen, the theory 
is unnatural. 

This definition assigns parameters as being natural or unnatural, but does not 
give a quantitative measure of the degree of naturalness. The current use of the term 
naturalness in the particle physics field has developed into referring to the sensitivity 
of observable quantities to small changes of the fundamental parameters governing 
the respective physics. 

To illustrate this identification of naturalness, a crude example would be walking 
into a room and observing a pencil standing on its point on a table. This situation 
requires fine tuning in its positioning (although much less fine tuning than present 
in the SM hierarchy problem) and so is considered to be in an unnatural state. It 
can be commented that we are able to control the positioning of the pencil, but we 
cannot change the Lagrangian parameters of a physical theory. However, if these 
parameters are actually governed by a more fundamental theory, then the analogy 
may be appropriate. 

The sensitivity definition of naturalness characterises the presence of apparent 
theoretical accidents required for observed physics. The degree of naturalness is then 
associated with the magnitude of the sensitivity. It is common to consider tuning of 
more than 1 part in 100 as unnatural, but being an aesthetic criterion, the transition 
from a natural to an unnatural theory is not immediate. 
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A quantitative measure first introduced by Ellis et a/ [S], but commonly referred 
to as the Barbieri-Giudice measure [9] in the literature is defined by 



^ 



p dX 
X dp 



:i.86) 



where X is some observable quantity, and p is some input parameter. 

For the case that the observable is a linear combination of the parameters, this 
measure extracts the contribution from parameter p and normalises with respect to 
the magnitude of the observable. If there exist large cancellations in order to generate 
the observable, then A is large and the observable can roughly be considered to 
require a fine tuning of 1 part in A with respect to the associated parameter. If each 
parameter contributes a similar amount, then A ~ 1/n where n is the number of 
parameters. This is the natural limit. 

In general, observables are not linear combinations of fundamental parameters. 
The above measure has the useful features that it is dimensionless by construction 
and is invariant under rescaling of the observable or parameter. 

The above measure does not always agree with the qualitative expectations of 
naturalness though. For example, the mass of the proton is extremely sensitive under 
this measure to small changes in the strong gauge coupling constant at a high energy 
scale. However, dimensional transmutation is proposed to be a very natural way of 
generating a mass scale. Some care must then be taken in interpreting the results of 
such a measure. 

Anderson and Castafio suggested a modification to the naturalness measure [TU] . 
such that the ratio of A^ to some appropriate average (Ap) over the parameter space 
is considered. A suitable probability distribution must then be assigned across the 
parameter space for the averaging process. This prior is not well-defined and in- 
troduces further ambiguity. It has been common to restrict the averaging integral 
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over parameter space to the experimentally viable region. For any given point in 
parameter space, the output of this fine tuning measure then depends on the current 
experimental status, which makes comparison with previous or future calculations 
non-trivial. 

This Anderson-Castafio measure will find reduced values as the experimental lim- 
its are strengthened, or after discovery as more precise measurements are made. This 
measure is thus more a description of likelihood than absolute sensitivity to funda- 
mental parameters. As the latter interpretation has been associated with naturalness 
here, the Anderson-Castano modification is not used for the analysis. 

An alternative modification to the Barbieri-Giudice measure was proposed by 
Ciafaloni and Strumia [llj specifically for experimentally observed parameters. This 
involves weighting the Barbieri-Giudice measure by (ap/p), where a^ is the experi- 
mental uncertainty in its value. The motivation is for identification of a probability 
with the inverse of the measure that the observed physics is described by the set 
of input parameters (within experimental errors). The weighting rescales the mea- 
sure to account for possible variations in the measured input parameter. This fine 
tuning measure shares some of the problems of the Anderson-Castano measure in 
interpretation of its results, and comparison with other work. However, when the 
parameters being tested involve a mixture of measured and unmeasured parameters, 
more importance is then attributed to the unverified part of the model. 

Beyond the choice of quantitative measure, there are also ambiguities in how to 
define the overall fine tuning of a given model. The option adopted for the following 
analysis is to take the maximum A^ over all choices for the parameter, p. This choice 
returns the worst sensitivity with respect to any parameter. Another option com- 



monly employed is to consider \/Yl ^p- Returning to the example of an observable 
that is a linear combination of underlying parameters, if only two terms involve a 
large cancellation, then the difference from the previous option is a factor of \/2. 
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These two choices of a maximum, or power sum of the individual A^, do not 
distinguish between few or many parameters requiring fine tuning. One may also 
consider the product of A to test for such cases. However, as the lower limit of A 
is zero and not unity the results of such a measure can be misleading. This section 
mentioned various limitations with all discussed fine tuning measures, and none are 
perfect for an indication of naturalness. An awareness of these limitations must be 
maintained when analysing the results of the following chapters. 
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Chapter 2 



Fine tuning of electroweak scale 
within minimal supergravity 



The discussion in the introductory chapter provides phenomenological and theoretical 
motivations for extending the Standard Model (SM) symmetry to include minimal su- 
persymmetry (SUSY). The Minimal Supersymmetric Standard Model (MSSM) must 
have some mechanism to spontaneously break SUSY in order to agree with the ob- 
servations of our universe. However, as argued in the introduction, phenomenological 
constraints force SUSY breaking to occur in some hidden sector. The communica- 
tion of SUSY breaking is here assumed to be mediated via minimal supergravity 
(mSUGRA) interactions. 

One of the main motivations for low energy SUSY is its use as a solution to the 
hierarchy problem. Having used naturalness to motivate the introduction of SUSY, it 
then must be checked that natural theories are experimentally viable. The measure 
can also be used to identify regions of parameter space that have an acceptable 
degree of fine tuning, and so determine promising search strategies for discovery of 
such physics. If naturalness is abandoned, then any new physics not excluded by 
experiment may be equally likely to be realised, which leads to immense freedom in 
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model building. 



Section 12m presents the physics of the Higgs potential and determines the one and 
two loop leading log corrections to the Higgs potential. Section fI72\ then evaluates an 
analytic formula for the fine tuning of the electroweak scale based on the former Higgs 



potential. Section [2^31 presents results for fine tuning across the mSUGRA parameter 
space. The dependence on constraints, input parameters and observables is discussed, 
and the influence of the fine tuning with respect to individual parameters analysed. 
The natural regions of the parameter space are then identified and naturalness limits 
obtained for the parameter space and spectrum. Finally the phenomenology of the 
natural realisations are discussed and a summary given. 

2.1 Higgs potential 

The introductory chapter identified the physics of the Higgs sector as contributing 
significantly to the degree of naturalness present in a model. It was also mentioned 
that this sector is soon to be tested by experiment and so is of topical interest. 
This section discusses the general structure for self interactions of the MSSM Higgs 
scalar fields at the electroweak scale, and will be referred to in following sections for 
explaining the behaviour of the naturalness measure. 

Introducing an alternative notation to that presented in Section [L6| with H^ = H^ 
and if 2 — H^i the scalar potential involving Higgs field terms is parameterised by the 
following, 

V = ml\Hi\'^ + ml\H2\'^-iTnlHi-H2 + h.c.) 

+ ijAiliiil +7^\2\H2\ + \3\H1\ \H2\ + A4 l-fii ■ -^21 



+ 



i As (i/i ■ H2f + Ae \H,\^ {H^ ■ H2) + A7 I//2I' (i^i ■ H2) + h.c. 



(2.1) 
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as for a general two Higgs doublet model, and H1.H2 = H^ H2 — H^ H2 ■ 

For analysis of the scalar potential, it is convenient to introduce the following 
parameters, 

rn? = ml cos^ /3 + m^ sin^ f3 — ml sin 2/3 
A = ^ cos^ /3 + ^ sin^ /3 + ^ sin^ 2/3 + sin 2/3 (Ag cos^ /3 + A7 sin^ /3) (2.2) 

where A345 = A3 + A4 + A5. On the Higgs field acquiring vevs, {Hf) = vsm(3/\/2, 
and {H2) = f cos/3/v2, a similar form to the single Higgs potential is found. The 
vacuum minimisation conditions provide the following constraints: 

v' = -^, 2A^ = m^f| (2.3) 



or equivalently. 



2m? 2 

— ^ — m, 



sin2/3 - .-n+"^2 + T [^1 4 + ^2 4 + ^345 + (\ + ^7)^2/3 + \ cot /3 + A7 tan /3] 

ml - ml tan^ /3 = - ^ [cos^ /3 ( A^ - A2 tan^ /3) + sin 2/3 (Ag - A7 tan^ /3) ] (2.4) 

In the two Higgs doublet model, there are eight real degrees of freedom. Elec- 
troweak symmetry breaking will adopt three of these degrees of freedom in order to 
give mass to the electroweak gauge bosons. This leaves two CP even fields, /i", H^, a 
CP odd field, A^ and two charged fields, H^. Considering the case of real parame- 
ters, the CP eigenstate fields do not mix. The masses of the CP even fields are given 
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below, using the notation of 



771, 



^-i (2A5 + X, cot f3 + X, tan f3) (2.5) 



''^h,H — 2 



m\ + v^ (2A + A5) ± V K + ^' (A5 - 2A)]' + Av^ Aj (2.6) 

^5 = ^ (-^1 + -^2 ~ 2A345) + A5 — ^ (Ag - A7) 

^6 = ^ (-^3451 C/3 - A3452 S^) + ^ (Ag + A7) + -5^ (Aq - A7) 

where A345J = A345 — A^, (s/3, c^) = (sin/5, cos/9), and /i° is conventionally lighter than 
H^. For large splitting (large m^), ml reaches an upper limit of 2Af^ (which tends 
to Ajf^ for large tan/3). In order to allow a large h^ mass, a large effective quartic 
coupling is then required. The masses of the A^, H^, H^ are all similar in this limit, 
with the h^ remaining relatively light. A decoupling limit is thus obtained where the 
light Higgs is then SM-like. 

For the MSSM, the only tree level contributions to the quartic couplings come from 
the electroweak D-terms. The effective quartic coupling is then, A = {g^/8) cos^ 2(3, 
and the individual A^ are given by: 

Ai = A2 = H^? + ^2), h = -A9l~9l), K = -\gl A5,6,7 = (2.7) 

where g^ ((72) is the hypercharge (weak) gauge coupling and g^ = gf + g^. 

Given that m\ = g^ f ^/4, the upper limit to the MSSM light Higgs is ?7i| cos^ 2/3 
at tree level. For SM-like couplings (found in the large m^ limit), this scenario is 
experimentally ruled out. It is then necessary to rely on the loop corrections in 
order to possibly escape this constraint through the increase of the effective quartic 
coupling. 



48 



2.1.1 1 Loop Leading Log (ILLL) Terms 

Using the MSSM renormalisation group equations (H], but neglecting the mass scales 
of the SM and superpartner states, the scalar potential remains of the form, 

\/(°) = ml\H,f + ml\H,f-mliH,H, + h.c.) + sl{\H,f-\H,fY (2.8) 

However, once low energy scales are reached, the masses of the MSSM states signifi- 
cantly affect the running. This can be approximately accounted for by running down 
to the mass scale of a given state, integrating that state out of the Lagrangian, and 
then running below the mass scale using the updated renormalisation group equa- 
tions. The effect of including the threshold dependence introduces corrections to 
the scalar potential which are accounted for by addition of the Coleman- Weinberg 
potential [13]. 

The 1 loop Coleman-Weinberg potential has the form, 

V^'^ = 6^E(-1)''^(2^^ + 1)^^^U^°S3I-|) (2.9) 

k 

where rrif, is the field dependent mass, defined as the second derivative of the potential 
with respect to the given field, evaluated for the vacuum state. The degeneracy factor, 
gi^, is 6 for (s) quarks, and Jf^ is the particle spin. All parameters in eq fl2.9j) are 
evaluated at the scale Q using the RGEs which ignore the particle masses. The field 
dependent squark masses are (neglecting 0((?^) terms): 



m^^ ^ Ml + h',\H,\' + 3^{\H,\'-\H,\')Th,\AH,-f,H*,\ (2.10) 
ml^ ^ Ml + hl\H,\' + ^{\H,f-\H,\')TK\\H,-f,H;\ (2.11) 



and where TnQ^u,Di.^s) — ^s i^ assumed. Here, Mg is the soft SUSY breaking squark 
mass evaluated at the squark mass scale. 
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One can expand the non-linear (logarithmic) field dependence in V^^^ in inverse 
powers of 1/Ms in order to find the dominant threshold corrections, which come from 
the third generation squarks: 



K 



(1) 



m-rr^ 



t {ht m' + 2h1 Ml \H,\' + hi \A,H, - f^Hll') 



+ h 



4 l^t^h^f;^^ 



H. 



12M| 



+ ai[\H,\'- \H,\^) 2t h] \H^\' + 2M| (t - 1) + 



Ml 



(2.12) 



V, 



(1) 



167r^ 



t {hi \H,t + 2hl Ml \H,\' + hi \A,H^ - f,H;f) 



+ rib m. l^il 



Mj- 



12M'§ 



+ ^-^{\H,f-\H,f){2thl\H,f + 2Ml{t-l) 






(2.13) 



where 



t = log{Ml/Q' 



(2.14) 



When running below the EWSB scale, the inclusion of higher dimensional terms 
(threshold corrections) lead to a re-summation such that Mg is replaced by a mass 
scale related to the physical particle masses [H] . For the results found using analytic 
formulae presented in this thesis, the geometric mean of the stop masses is used in 
the place of Mg. 

The above equations are valid down to the top mass scale; below this scale thresh- 
old corrections from the top quark should also be included. The dominant effect of 
running below the top scale can be absorbed by setting Q in the above equations as 
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the "running" top mass evaluated at the scale Q instead of the pole mass. 

From eqs. (12. 8p to fl2.13p one obtains the parameters entering in the scalar potential 
( I2.ip . evaluated at the scale Q (below M5), in the one- loop leading log approximation: 
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These analytic results agree with [15] when ignoring the stop mixing corrections to 
the D-terms. 
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2.1.2 2 Loop Leading Log (2LLL) Terms 

The two-loop leading log (2LLL) Coleman- Weinberg potential can be found in the 
arXiv version of [I6] to 0{gi h^, gl hi) and 0(/i^, /i^ hi, hf hf, hi), see also [HI M\ for 
the general case. The method of the previous section may be used to determine the 
2LLL contributions to the Higgs scalar potential, however here a similar approach 
to that in [H] is used, to RG-improve the 1-loop result into a 2LLL result. A step 
approximation is applied to the /3-functions so that the MSSM RG eqs are used 
between the GUT and stop mass scale, then the 2HDM SM RG eqs between the stop 
and top mass scales, and finally the top is integrated out to reach the electroweak 
scale. 

When setting the renormalisation scale in eqs fl2.15p to f l2.23p as Q = Mg, the 
logarithmic terms are removed but the finite corrections from stop mixing remain. 
These results are then used as boundary conditions for the parameters at the scale 
Mg. A series expansion of the RG eqs can then be applied: 

XiQ) ^ XiMs)-PAMs)t + lP[iMs)t' + 0{t^) (2.24) 

= A (Ms) - /3a (g) t - i (3[ (Q) e + {e) (2.25) 

where (3^ = dp/dlogQ^. Eventually, all parameters will be expressed at a scale Q as 
in the Coleman- Weinberg potential approach. For a /3^-function of the form bX + c, 
eq (I2.25P becomes 

A ^ X{Ms)-t[bX{Ms)+c\ + t^[bc-lP[ + 0{X)] (2.26) 

where the couplings are evaluated at the scale Q unless stated otherwise. The (3- 
functions for the 2HDM SM [T^] are listed below, neglecting 0{h'^) terms, and with 
the /9^ -functions also neglecting 0{g^,g'^X^, Xf) terms: 
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and finally 






(2.29) 
(2.30) 



Using (12.261) . the analytic 2-loop results in fT5] are then recovered when the same 
level of approximation is considered. For example, 



3fe? 



[X,{Ms)-X,a^t]-h,t+ a^h, + ^{2P^.-P^]+0{X) t^ (2.31) 



6?i2 



3/if 



A2 {Ms)-X,^t + 1^ [t + ^ (I /^^ + i /.,2 - 8^1) e] (2.32) 



The couplings entering in the expression of Aj {Mg) are re-expressed in terms of 
their values at the scale Q, (with a logarithmic correction which compensates for the 
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running below M^): 



ht{Ms) 
ht{Ms) 



ht{l + ji^{9e,+hl-lQgi)+0{gH,t')) 
ht{l + ji^ {9hl + hl- IQgl) + 0{gH,e)) 



(2.33) 
(2.34) 



This leads to the following expression, in agreement with [15], when the stop mix- 
ing contributions to the D-terms in the potential are again neglected (which is not 
assumed for the analytic results presented in the following sections): 



A. 



^(l-S^)-S^^[l + I^(%^-5/^?-16^I)] 









^^ + ^ + I^(¥ + f-8^l)(^t + ^) 



(2.35) 



where the following notation has been introduced. 



^t. 






Mi 



12Mi 



(2.36) 



Note that these results assume that the CP odd Higgs mass is not decoupled. If this 
is the case, the usual SM /3-functions should be used. The effective quartic coupling 
at the EW scale when m^ < M^ is given by: 



A ^ ^cos^2l3 



l-Tt{hl + h', + {hl-h',)sec2f3)t 



+ 



3hf . 4 n 



3/1? , h'i 



t+^+Tt^r-f+Y-^ai w+n+6. 



I 3/i* 4 n 

+ I6^ COS /3 



t + ^ + Ti^m^ + '^-89l)(x,t + t^)+6., 



(2.37) 



with the following notation: 
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where X^^ is defined as X^^jy^A^y^ — )■ A^^) with 

A^ = A^ — n cot P 

if, = Afc-/itan/3 (2.40) 

A similar but distinct result is obtained when m^ ~ Mg (notably no 5^ terms and 
a different dependence on tan/3 and the mixed Yukawa couplings). The threshold 
corrections are dependent on where the CP odd Higgs decouples. The same procedure 
is applied to determine the 2LLL threshold corrections to the mass terms. The results 
using these analytic formulae for the loop corrections assume that the CP odd Higgs 
has not been decoupled above the scale Q = niz- This is inevitably wrong, but the 
errors involved are found to be insignificant given the other approximations made. 

2.2 Analytic formula for fine tuning 

The fine tuning measure chosen for this analysis is the Barbieri-Giudice measure 
introduced in Section [T75] 



A = max A , Ap 



9 In 1)2 



d\np 



(2.41) 



where the observable is taken as the electroweak vev squared. This is not actually an 
observable, but this approach is nearly equivalent to taking an electroweak gauge bo- 
son mass as the observable, where equivalency is broken by small quantum corrections 
to the self energy and gauge coupling running. It will turn out the approximations 
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applied for calculation of the fine tuning generate more uncertainty than these choices 
of "observable" introduces. 

The fundamental parameters, {p}, are considered to consist of the soft SUSY 
breaking parameters and the Higgs bilinear parameter within the superpotential, 
{rriQ, Tn\i2i v4q, i?Q, /Xq}. The dominant dimensionless parameters controlling the renor- 
malisation group equations, the strong gauge coupling and top Yukawa coupling, 
{03, h^}, were also tested for fine tuning. However, the latter results were not consid- 
ered as genuine measures of naturalness and were not included when finding max A 
. For example, the influence of the top Yukawa coupling leads to radiative symmetry 
breaking which, as will be discussed later, produces similar sensitivities as found for 
the dependence of the strong coupling on the proton mass scale. The Barbieri-Guidice 
measure is then not considered appropriate for these dimensionless parameters. 

As gravitational interactions mediate supersymmetry breaking in mSUGRA, the 
universal pattern of soft masses determined in Section 11.51 is valid around the Planck 
scale. These are then considered to be the fundamental parameters for the fine 
tuning tests. One motivation for universality relied on gauge unification. However, 
if unification occurs then the running of the Lagrangian parameters between the 
unification scale and Planck scale is certainly modified with respect to the defined 
MSSM case. The apparent scale of gauge unification in the MSSM is ~ 2 x 10^^ GeV. 

The running of soft masses between the unification and Planck scale will be ig- 
nored in the analysis. The relatively small hierarchy of these energy scales, and a 
preserved universality within any gauge unified multiplet may make this a reasonable 
assumption. The physics below the gauge unification scale is then described by the 
MSSM, and the associated renormalisation group equations are used to determine 
the low energy physics and their sensitivity to variations of the high energy inputs. 

Note the choice of testing a mass- dimension- squared observable to sensitivity in 
mass-dimension-squared parameters for the SUSY parameters. This may be consid- 
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ered appropriate, as in the tree level vacuum minimisation condition, f ^ = im? /\. The 
effective quartic coupling does not depend on the SUSY parameters at tree level, and 
so the observable is a linear combination of the chosen input parameters. As argued 
in Section ll.8[ the measure then extracts the sole term within this linear combination 
that depends on the tested input parameter, and normalises the result with respect 
to the observable's magnitude. If the mass dimensions selected were not equivalent, 
then the measure would be scaled by a constant factor. 

In general, the loop corrected mass parameters and quartic couplings in the two 
Higgs doublet potential will depend on {/3,p, f }. The measure is then given by. 



v^ dp 



1)2 



A2 \dv'^ 9/3 dp "^ dv^ dp 
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X \ dp dp dp dv'^ dp 



(2.42) 



where dP/dp can be determined by differentiating the second minimisation condition: 
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(2.44) 



Using the vacuum minimisation equations to replace terms, the following formula 
for the measure is obtained. 
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(2.45) 



This agrees with the analytic result by Casas et al [21] when ignoring the vev depen- 
dence of the Lagrangian parameters generated by loops (by setting r = 0). 

The impact of including the r factor tends to reduce the fine tuning by 0(30%), 
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however it can also increase the fine tuning for certain points in the parameter space. 
In the following analysis that is concerned with the minimum fine tuning obtained 
from scans of the parameter space, the effect of ignoring r simply scales the curves 
of minimum fine tuning found in the subsequent plots by this 0(30%) factor. The 
following analytic results are presented for the case of r = 0, which is equivalent to 
considering the superpartner mass scale fixed at which they need to be integrated 
out. The analytic results are not used here for determining high precision fine tuning 
values, but as an indication of the preferred regions of parameter space. The final 
results are determined numerically via a finite difference method within SOFTSUSY, 
and so do not use this analytic formula or its assumed approximations. 

In the case of ignoring the parameter dependence of the quartic couplings (valid 
at tree level when the input parameter is not an electroweak gauge coupling), the 
"master formula" of [22] (see also [231 [21]) is also found. 



A. 



(tan2/3-l)m| 



dp 



tan/3 
cos 2/3 



tan^ /3 ^ 



1 + 



nir, 



m^+m 



op ' \ op op 



(2.46) 



This formula is sometimes used as the starting point in analyses that evaluate elec- 
troweak scale fine-tuning, by using in it the loop corrected soft masses. However, for 
accurate estimates, it is necessary to take full account of radiative corrections. 
In the limit of large tan/3, eq (I2.45P reduces to 



A _ [4A7 {mlY-i\'^7nl] + [\'^v-^+2{mly] [A3452+2K-mi)/r,2] x.^/ , n^ 



2Af d2-A2 [A3452 t>2+2(mf-m: 



(2.47) 



—)■ 



\^v^ 



2 {ml)' + A'2 v^ + Av^ / "^^fJlT^t^S, 

\ 2.J ' 2. \^A3452 i>2+2{mf-m.^) 



(if IA7I < |A2|,|A3452|) 



where x' = dx/dlnp is the partial derivative of x wrt In p. 

In order to match physical observations, the vacuum minimisation conditions are 



used to set the magnitude of /i, using the electroweak gauge boson mass as input. The 
sign of /i however remains a free parameter. It is also common to exchange one of the 
mSUGRA parameters, Bq, by tan/3 in order to characterise the position in parameter 
space. Although, as tan /3 and m^ are not fundamental parameters of the high energy 
Lagrangian, the fine tuning will still be evaluated with respect to variations in Bq 
and fi. 

2.3 Fine tuning results from parameter space scan 

The following presents results of evaluating the fine tuning measure over an experi- 
mentally allowed region. The fine tuning was first always evaluated using the analytic 
formula of the previous section, including the dominant third generation supersym- 
metric threshold effects in the scalar potential. The evaluation of the analytic formula 
was completed using Mathematica with a routine optimised to run quickly. The re- 
sults were then used to identify regions of parameter space that possessed small fine 
tuning. 

For the cases where an analysis of the general behaviour of fine tuning was con- 
sidered sufficient, the analytic results were then used. This included testing the 
importance of various approximations. For the cases where more precise results were 
desired, the fine tuning measure was re-evaluated using the numerical routines of 
SOFTSUSY 3.0.10 ^25j. It proved computationally prohibitive to begin a search of 
the parameter space using the numerical routines. This approach would also have 
made analysis of the influencing factors much more challenging. 

The mSUGRA parameter space scanned over was selected to satisfy the following 
constraints: 

• radiative electroweak symmetry breaking (EWSB) is found 

• non-tachyonic SUSY particle masses 
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(to forbid colour and charge breaking (CCB) vacua) 

• experimental constraints satisfied 

• if stated in plots, consistency of rrih with the LEPII SM-like bound (114.4 GeV) 
where the experimental constraints include the bounds on superpartner masses, elec- 
troweak precision data, b physics constraints, and the anomalous magnetic moment 
of the muon, as detailed in Table 12.11 



Constraint 


Reference 


SUSY particle masses 

6a^ < 366 X 10"" 

3.20 < 10^ Br(6 -^ 57) < 3.84 

Br(5, -^ fifi) < 1.8 X 10-8 

-0.0007 <5p< 0.0012 


Routine in Mi cr OMEGAS 2.2, "MSSM/masslim.c" 

PDG (sys. and stat. Icr errors added linearly) 

PDG (sys. and stat. la errors added linearly) 

Particle Data Group (PDG) 

Particle Data Group (PDG) 



Table 2.1: Constraints tested using MicrOMEGAs 2.2. PDG: http://pdg.lbl.gov/. 

Regarding the Higgs mass constraint, the values determined for a given point 
in parameter space with SOFTSUSY agree with that found using SuSpect [26] within 
0.1 GeV, but can differ by 2 GeV [27] from the value found using FeynHiggs [28j. The 
difference is a 3-loop effect caused by different renormalisation schemes. The value 
of the Higgs mass obtained from SOFTSUSY should then be considered to possess a 
theoretical uncertainty of 0(2 GeV). For the numerical results, the SOFTSUSY Higgs 
mass is used. 

Fig 12.11 presents results from the analytic fine tuning measure, demonstrating 
generic features for the minimum fine tuning allowed within the parameter space 
region scanned, which are discussed in the following subsections. 



2.3.1 Dependence on Higgs mass 

For any fixed tan/3. Fig 12.11 shows similar behaviour in the minimum fine tuning 
curves vs Higgs mass. There exists a global minimum at an intermediate Higgs mass, 
with rapidly increasing fine tuning as the Higgs mass is increased or decreased. It 
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Figure 2.1: Minimum (analytic) fine tuning versus 2LLL Higgs mass, showing the influence 
of loop effects and tan/3 on fine tuning. All constraints listed in Table [2TT] are included. The 
upper and lower lines associated with the coloured regions are the 1-loop without thresholds 
for A and soft masses and "full" 1-loop results respectively (similarly for the grey region, 
for all tan/3). The minimum 2-loop fine tuning is found between these two cases. The solid 
lines refer to the scan 2 < tan /? < 55. The black points give the positions of minimum A 
for fixed tan/3 from 2 to 4 inclusive in steps of 0.5. 



has been previously argued that if SUSY is a solution to the hierarchy problem then 
the superpartners that affect the Higgs potential (eg stops) are preferred to be light. 
The non-observation of superpartner states has introduced some tension by requiring 
heavy masses. The least fine tuning is then expected when the stops have masses 
just above the experimental limits. This corresponds to a certain value for the Higgs 
mass, for a given tan/3, which can be read from Fig |2.1[ 

In order to get lighter Higgs mass, the A^ ^ or /i parameters need to be large with 
an appropriate sense in order to reduce the loop contributions. However as the A 
measure grows like A^ or /i^ (at least before leading log corrections, see Section [2. 3. 6p . 
this leads to increased fine tuning. 

On increasing the stop masses above their experimental limit, the Higgs mass is 
increased. Inverting the Higgs mass formula near the decoupling limit gives, m^ oc 
exp(m^). As the stop-masses-squared are proportional to the soft SUSY breaking 
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mass-squared terms at tree level, the A measure is then expected to roughly scale 
exponentially as the Higgs mass is increased. 

2.3.2 Dependence on tan/3 



Fig 12.11 also shows a clear dependence of minimum fine tuning on tan/3. As tan/3 is 
increased, the tree level Higgs mass is increased through an increase in the effective 
quartic coupling. The analytic formula for fine tuning exhibits an inverse proportion- 
ality to this effective quartic coupling. The global minimum fine tuning for a given 
tan /3 is then expected to scale as m^^ evaluated at this minimum, and this is roughly 
observed. 

2.3.3 Dependence on N-loop leading log corrections 

Similarly for the dependence on tan/3, the infiuence of the leading log corrections to 
the effective quartic coupling shifts both the Higgs mass and fine tuning. The ILLL 
terms tend to increase the quartic coupling by a factor of 0(2), reducing the fine 
tuning by a similar factor and increasing the Higgs mass by the square root of the 
proportional change. This behaviour agrees with the results of [32] where a limited 
parameter space scan was completed for the experimental constraints as of 1998. 

The inclusion of 2LLL corrections is seen to increase the fine tuning with respect to 
the ILLL results. This again follows as the loop diagrams with strong gauge coupling 
dominate, and these lead to a reduction in the effective quartic coupling. The 2LLL 
corrections similarly decrease the Higgs mass. 

As the change in fine tuning from 2LLL effects can introduce 0(2) corrections for 
a fixed Higgs mass, they are significant in evaluating naturalness of the parameter 
space. The 2LLL results are now used throughout the following analysis, and when 
numerically determining the fine tuning via SOFTSUSY, the 2-loop renormalisation 
group equations will be used. 
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The effect of going to 3-loop calculations generally introduces 0(1%) changes to 
the Higgs mass. Having observed a correlation between minimum fine tuning and 
Higgs mass, it is then reasonable to expect that the fine tuning may vary by 0(2%) 
at this level. This error is well within that generated by the various approximations, 
and so it is not currently worthwhile to include 3-loop corrections for evaluation of 
the fine tuning. The interpretation of such small differences in fine tuning is also not 
so meaningful. 

2.3.4 Dependence on superpartner mass limits 

Fig 12.21 demonstrates the effects of varying the superpartner mass limits. This is 
useful to compare with previous fine tuning results when experimental constraints 
were weaker, and understand how stronger mass limits can affect the minimum fine 
tuning allowed. 

It is found that with the current experimental limits, the light chargino mass con- 
straint is almost wholly responsible for restriction on the minimum fine tuning found 
across the parameter space due to superpartner mass limits. This is a consequence of 
recent experiments, as when older experimental limits are used, the neutralino mass 
limit is found to be dominant. In both cases, it is masses controlled by the gaugino 
soft mass term, 171-^^,2 that are found to be most significant for limiting fine tuning. 
The gluino mass limit is close to providing a matching constraint on fine tuning as 
the chargino limit, but is currently sub-dominant. 

The current chargino mass limit varies according to the position of parameter 
space being tested, but for many cases, m ± > 94 GeV is roughly imposed. On 
decreasing the chargino mass limit, the minimum fine tuning is seen to reduce in 
Fig 12.21 for m^ < 114 GeV, with little effect for heavier Higgs masses. On increasing 
the superpartner mass limits controlled by Tn-^,^, there is a roughly parallel movement 
in the fine tuning vs Higgs mass space of the minimum fine tuning curve at low Higgs 
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Figure 2.2: Minimum (analytic) fine tuning vs Higgs mass, showing the influence of mass 
constraints on fine tuning. The results are at 2-loop with the upper shaded (coloured) areas 
connecting the case of only applying the SUSY spectrum constraints (lower line) to that 
with all constraints listed in Table 12.11 (upper line) . The lower shaded (coloured) areas 
connect the cases of only applying a chargino lower mass limit of 80 and 94 GeV for the 
lower and upper lines respectively. The results for the scan 2 < tan/3 < 55 are also shown 
by the grey shaded area, with similar convention for upper/lower continuous lines delimiting 
it. 



mass. At large Higgs masses, the exponential increase in fine tuning is still governed 
by large stop masses well above the experimental limits. 

Note that this can be used as a test for naturalness. If a limit for fine tuning of 
1 part in 100 (A = 100) is demanded, then upper mass limits on superpartner states 
and Higgs mass can be obtained. If any state is observed in breach of these limits, 
then it can be concluded that a "natural" mSUGRA theory cannot be realised. This 
may point to requiring an alternative theory to mSUGRA, or the abandonment of 
the principle of naturalness. The idea of setting of naturalness limits will be returned 
to when discussing the numerical fine tuning results. 
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Figure 2.3: Minimum (analytic) fine tuning vs Higgs mass, showing the influence of the h — )• 
S7 constraint. The results are at 2-loop with the lighter shading connecting the case of only 
applying the SUSY spectrum constraints (lower line) to that with also the 6 — )• 57 constraint 
listed in Table 12.11 (upper line) . This upper edge of this shading is indistinguishable from 
the solid line which includes all constraints in Table 12.11 The darker shading extends up to 
the minimum fine tuning limits for the stronger case, 3.52 < 10^ Br(6 — )• 57) < 3.77, with 
the other constraints as given in Table 12. 1[ The results for the scan 2 < tan /? < 55 are also 
shown, between the two continuous and almost parallel lower curves. 

2.3.5 Dependence on 6 — > 57 constraint 

Figure 12.31 gives the impact of the 6 — )■ 57 constraint on A. The lower limit of the 
6 — )■ 57 constraint for a given coloured area (fixed tan/3) restricts the right hand 
edges of the plot, while the upper limit restricts its left hand side. These curves also 
depend on the mass limits - these constraints are not fully independent. 

For the experimentally allowed area of mh > 114.4 GeV, the impact of the b ^ s'-f 
constraint is rather small. The combination of the SUSY mass limits and 6 — )■ 57 
constraint currently dominate the restriction on how small the fine tuning could be. 



The other constraints listed in Table 12.11 do eliminate further mSUGRA points, 
but have a negligible effect on the fine tuning limits. With the current mass limits, a 
change in the 6a constraint by factors of 2 or more does not noticeably affect these 
results. 
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Relaxing the 6 — ?■ 57 constraint to be within the 3a experimental error, or even 
neglecting the constraint totally only introduces 0(10%) changes in the minimum fine 
tuning allowed for the full parameter space scan. This was verified with the numerical 
results also. However, for fixed large tan/3, the importance of this constraint remains 
significant. 

2.3.6 Dependence on the SUSY Lagrangian parameters 

The remaining analysis now includes a re-evaluation of the fine-tuning using SOFTSUSY 
in the regions of parameter space that were found to populate the low fine tuning 
points in each plot. The SOFTSUSY results were verified to exhibit the same generic 
features, with the quantitative differences attributed to a more robust evaluation of 
the renormalisation group equation solutions. 

The results of this section identify which regions of parameter space permit small 
fine tuning of the electroweak scale, and will later be used to comment on the phe- 
nomenology of natural mSUGRA theories. In Figure |2^ the dependence of the total 
fine tuning with respect to the Lagrangian parameters is plotted. 

To understand the behaviour of the electroweak fine tuning, it is useful to consider 
one of the vacuum minimisation equations, which when ignoring quantum corrections 
to the quartic couplings simplifies to: 

2 ~ tan2/3-l l^-^^J 

Using 2-loop RGE solutions which neglect MSSM state masses for tan/3 = 10, 

ml{mz) ^ 0.99 ;u^ + 0.946 mg + 0.331 m^/2 + 0.044 Aomi/2- 0.013 A^ (2.49) 
ml{mz) ^ 0.99 ^^ - 0.080 m^- 2.865 m^2 + 0.445 Aomw2- 0.099/1^ (2.50) 
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(d) Fine Tuning vs ttIq 
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Figure 2.4: Dependence of minimum fine tuning on SUSY parameters (// > 0). The solid, 
dashed and dotted hnes are as explained in Fig I2.5[ No bound on m/j is applied in figure 
(a). In (c), (d), (e), the darker shaded regions are eliminated when ?7i^ > 114.4 GeV is 
applied for the case with the central {a^^ra^) values. In (b) and (f), ra^ > 114.4 GeV is 
applied, and the points in (f) are only for the central {a^,m^) values. 

The evaluation of A 2 for example, is roughly equal to (2/iQ/?7i|), when neglecting the 
leading log corrections. Small values for the Lagrangian parameters are thus preferred 
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to obtain small fine tuning. The quadratic dependence on the respective Lagrangian 
parameter is a common feature, which is observed in the plots for large parameter 
values. 

It is the large cancellation between the /Iq and fri^i^ terms that is often responsible 
for the large fine tuning (note however that this argument ignores the impact of 
quantum corrections to quartic couplings, known to reduce the fine-tuning). This 
leads to the approximate relation A 2 ~ A 2 • The rise in fine tuning at small mi/2 
is a consequence of the constraints such as the chargino mass limit. 

The near flat distribution of minimum fine tuning in mg is a result of the coefficient 
of ttlq in 1712 being driven close to zero. The fine tuning with respect to itlq then rarely 
dominates, until we reach values of rrih above the LEPII bound {rriQ at the edge of 
focus point region). The result of applying the Higgs mass constraint also excludes 
a region with small mi/2 at ttiq below 1.5 TeV. The focus point at tuq ~ 1.5 TeV 
where the minimum of mi/2 is possible, corresponds to the point where fine tuning is 
minimised. This only occurs for large tan/3, and this available "dip" in fine tuning 
in the mSUGRA space disappears as tan/3 is reduced. 

Figure [231(c) indicates that a small trilinear coupling \Aq\ < 1 TeV is preferred 
for the smallest fine tuning. This follows from a similar argument for preferring small 
mi/2- Increasing I^qI requires larger cancellations with /i to set the electroweak scale. 
However, once the Higgs mass constraint is applied, A^ is driven negative for small 
tan /3 in order to maximise the stop mixing. The related increase in the minimum fine 
tuning from being in this region of parameter space then follows. This is important 
for small tan/3 where the tree level Higgs mass is smallest. The sign structure of the 
UV parameter coefficients in 7712 leads to a preference in a small, positive A^. 
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100 no 

Higgs Mass /GeV 

Figure 2.5: Fine tuning vs Higgs mass, in a two-loop analysis. The data points are for 
2 < tan/3 < 55. The solid line is the minimum fine tuning with central values (03,771^) = 
(0.1176, 173.1 GeV). The dashed line corresponds to (03, mj) = (0.1156, 174.4 GeV) and the 
dotted line to (0.1196, 171.8 GeV), to account for la experimental errors in 03 and top mass 
|33j . This is the "worst" case scenario, when such deviations combine such as to give the 
largest variation of A. The LEPII bound of 114.4 GeV is indicated by a vertical line. Note 
the steep (~ exponential) increase of A on both sides of its minimum value situated near 
the LEPII bound. 

2.3.7 Dependence on strong gauge coupling and top mass 

The large quantum corrections to the Higgs mass provides a sensitive dependence on 
the supersymmetric physics, and therefore fine tuning. We now return to examination 
of fine tuning in the Higgs mass space. Fig 12.51 demonstrates how the fine tuning of 
our observed physics varies according to the experimental input for the strong gauge 
coupling and top mass at the electroweak scale. Their parameters were varied within 
their la experimental bounds to determine the combination that most strongly affects 
the fine tuning results. 

Note that this variation is philosophically different to testing sensitivity of ob- 
served physics to changes in the fundamental parameters. The variation found is not 
directly an indication of fine tuning, but a test of how accurate the fine tuning results 
can be considered to be. 

Varying the gauge coupling and top mass influences the leading log corrections 
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to the Lagrangian parameters. These affect both the Higgs mass and fine tuning 
to a similar order of magnitude. The experimental constraints are also affected, as 
the SUSY mass spectrum is generally modified. Also, the region of parameter space 
where radiative electroweak symmetry breaking occurs depends on these parameters. 
An increase oiaslmz) or reduction o{int{mz) by la have similar effects, which can 
be also understood from the relation between the mass of top evaluated at niz and at 
rrit. Keeping either a^ or mt fixed to its central value and varying the other within la 
brings a curve situated half-way between the continuous line and the corresponding 
dashed or dotted line. 

2.3.8 Individual fine tunings 

The individual contributions to A are shown in Figure 12.61 Below the LEPII bound, 
detailed calculations show that the minimal value of A is dominated by A^2 and 
this increases rapidly for decreasing rrih- For values of uih above the LEPII bound, 
A is dominated by ^^n?. This happens at the edge of the focus point region. The 
transition from the dominant A^2 regime to the dominant A^2 regime occurs near 
the LEPII bound value, where the fine tuning happens to be minimised. This is a 
consequence of the present day experimental constraints. 

The large values of A^2 (using the Barbieri-Guidice measure) when m^^ > 105 GeV 
suggests an apparent large sensitivity to the top Yukawa coupling in the region argued 
to have small fine tuning in the other plots. The small fine tuning with respect to 
all other parameters is a consequence of being near the scalar focus point, where the 
sensitivity to the scalar soft mass parameter is reduced due to the renormalisation 
group equations. 

To illustrate the presence of the scalar focus point, consider a moderate tan/3 
scenario with Aq = 0. The running of the soft mass terms is dominated by the top 
Yukawa coupling, where neglecting corrections from gauge interactions and subdom- 
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Figure 2.6: The plot displays the individual contributions A , to the minimum electroweak 
fine-tuning A presented in Figure 12.51 The largest of these for all nih gives the curve 
presented in Figure [231 At low rrih, A 2 (red) is dominant, while at large nih-, A^2 (blue) 
is dominant (with A^2 reaching similar values near 120 GeV). The transition between the 
two regions is happening at about 114.5 GeV. Note that in this plot the LEPII bound is 
not imposed at any time. Although A^2 (purple) is presented above for illustration, this 
contribution is always sub-dominant if using the Ciafaloni-Strumia definition of fine-tuning 
|llj . for measured parameters. In all other plots, A^2 is ignored for the calculation of A as 
it is not considered a genuine measure of naturalness, argued in the text. 



inant Yukawa interactions, the 1-loop RGEs are: 
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which leads to the following solution for the up-type Higgs soft mass term, 



m]iAh) 



i[3. 



l + iXj 



X 



R 



)(e^-l)]mi(to) 



(2.52) 



where m]^ J^t^) = {l + x^j^) m]j^{Q, P = //^^ g| dt, and t^ = log (Qi/Qo) for energy 
scale Qj. If the MSSM RGEs continued to arbitrarily low energy scales, there would 
always be a scale at which the up-type Higgs soft mass would be (approximately) 
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independent of the high energy scale input. However, observables must be indepen- 
dent of renormahsation scale. In this discussion, the mention of a scalar focus point 
implicitly implies that this focus point scale occurs at the electroweak scale where the 
vacuum state is set. For this scenario, the electroweak physics is then insensitive to 
the fundamental scalar soft mass. 

To avoid phenomenological problems in the MSSM, the high energy soft masses 
cannot be tachyonic. The top Yukawa coupling drives mfj to decreasing values as 
the energy scale is reduced below the gauge unification scale. Electroweak symmetry 
is then broken spontaneously by radiative corrections [31] as rajj becomes negative, 
perhaps requiring contributions neglected in eq f l2.52p . It is the same radiative cor- 
rections that then control the positioning of the focus point scale. 

In the mSUGRA scenario, a^^ r = 0, and the scalar focus point for low to moderate 
tan/3 is found where e^ ~ 1/3. Higher loop corrections, thresholds and neglected 
contributions modify this relation, but the observed top mass generally places the 
focus point scale close to the electroweak scale with better matching at large tan/3. 

Note that universality is not required for a focus point, as shifts of x^^^ — )■ 
x^ ^ + e^ ^ where e^ = — e^ will maintain independence of the high energy value. 
The ratios must still be fixed with one underlying mass-scale parameter and not co- 
incidental. Independent shifts in the ratios of other high energy scalar soft masses 
may also be found in the large tan/3 cases where the other third generation Yukawa 
couplings become significant [35]. However, for alternative mechanisms of mediat- 
ing supersymmetry breaking such as gauge or anomaly mediation, the pattern of soft 
masses generally does not conform to that required for a focus point at the electroweak 
scale. 

Returning to discuss the large apparent fine tuning to the top Yukawa coupling in 



Fig 12. 6[ one has to be careful whether this large result is meaningful. In the example 
of the QCD scale, there was a large value for A with respect to the high energy 



72 



gauge coupling input due to an exponential dependence for setting the confinement 
scale generated through running. Similarly here, there is an exponential dependence 
on a dimensionless parameter which controls the electroweak scale through radiative 
symmetry breaking. 

For the following analysis, radiative EWSB is considered a natural mechanism. 
The large result of the A measure with respect to the top Yukawa coupling is then 
attributed to an inappropriate choice of fine tuning measure given how this param- 
eter infiuences the electroweak scale, and ignored. Anyway, once electroweak fine 
tuning is at a level of 1 part in 100, there are also problems such as the hierarchy 
of Yukawa couplings that should be addressed. Further new physics may provide 
reduced sensitivities to the true fundamental parameters. 

2.3.9 Naturalness limits 

As mentioned earlier, the fine tuning measure A can be used to establish the remaining 
parameter space of mSUGRA compatible with a solution to the hierarchy problem. 
Assuming that A = 100 is the upper limit beyond which we consider that SUSY 
failed to solve the hierarchy problem, the following naturalness bounds are obtained: 

mh < 121 GeV 5.5 < tan/3 < 55 

/i < 680 GeV 120 GeV < mi/2 < 720 GeV (2.53) 

mo < 3.2 TeV -2.0 TeV < Aq < 2.5 TeV 

Note that the uncertainty in top mass and strong gauge coupling produces some 
variation in these limits. The Higgs mass limit varies at the 0(1%) level, and the 
AqjIti-^^,^, fi limits at the 0(10%) level. The mp limit is particularly uncertain though 
due to movement of the boundary where electroweak symmetry breaking occurs in the 
parameter space. The mg limit ranges from 2.6 to 4.2 TeV on varying the Standard 
Model inputs within their la uncertainties. 
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The naturalness restriction on the mSUGRA parameter space also leads to a 
restriction on the MSSM superpartner spectrum. The results so far demonstrate that 
electroweak fine-tuning has a strong sensitivity to parameters such as fi, mi/2, with 
a preference for lower values. Regarding the ttiq dependence, A has a rather fiat 
dependence. The states that are dominantly controlled by the yU, mi/2 parameters 
are then the most important in determining the naturalness of the proposed theory. 
These include the neutralinos, charginos and the gluino states. If any of these states 
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Table 2.2: Upper mass limits on superpartners in GeV such that A < 100 remains possible. 
These limits apply individually on the respective state, independent of other states. 



have masses in excess of those given in Table 12. 2[ it will require less than 1% tuning 



(A > 100) for the MSSM. These upper mass limits scale approximately as ^yK~, 
so they may be adapted depending on how much fine tuning the reader is willing to 
accept. 

In the limit of optimum naturalness, the fine tuning of the electroweak scale is 
found to be minimised at A = 8.8, for a Higgs mass of 114 ± 2 GeV. The uncertainty 
in Higgs mass is generated from the theoretical approximations, and the fine tuning 
may be considered accurate to 0(30%) with the given definition. These results are 
calculated at the 2 loop level. As noted previously, the most natural Higgs mass for 
mSUGRA theories is fixed by the current bounds on the superpartner spectrum, and 
not directly the LEPII Higgs mass bound. 

One may interpret the SUSY parameters corresponding to the smallest fine tuned 
region as being the most likely, given our present knowledge, and so it is of interest 
to compute the SUSY spectrum for natural parameter choices as a benchmark for 



future searches. This is presented in Table 12.31 where it may be seen that there are 
very heavy squarks and sleptons and lighter neutralinos, charginos and gluinos. This 
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has similarities to the SPS2 scenario 
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Table 2.3: A favoured mSUGRA spectrum with A = 15. Masses are given in GeV. 



2.3.10 Predictions for SUSY searches at the LHC 

It is clear that there is still a wide range of parameters that needs to be explored 
when testing mSUGRA. Will the LHC be able to cover the whole range? To answer 
this, note that for a fine tuning measure A < 100, one must be able to exclude the 



upper limits of the mass parameters appearing in Table 12.21 Of course the state that 



affects fine tuning most is the Higgs scalar and one may see from Figure 12.51 that 
establishing the bound rri/j > 121 GeV will imply that A > 100. However the least 
fine tuned region corresponds to the lightest Higgs consistent with the LEPII bound 
and this is the region where the LHC searches rely on the /i — ?■ 77 channel which 
has a small cross section and will require some 30 fb^^ at ^/s = 14 TeV to explore. 
Given this, it is of interest to consider to what extent the direct SUSY searches will 
probe the low fine tuned regions. Following the discussion in the previous section, 
the most significant processes at the LHC will be those looking for gluinos, winos and 
neutralinos. 

Studies of SUSY at the LHC [37] have shown that the LHC experiments have a 
sensitivity to gluinos of mass up to 1.9 TeV for ^/s = 10 TeV, and up to 2.4 TeV for 
^/s = 14 TeV and luminosity 10 fb^^. This corresponds to probing up to A = 120, 180. 
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In a previous study of fine tuning |2B] for fixed Aq = 0, tan f3 = 10, M^ = 174 GeV, 
a scan over the {171^,171-^^,^) space identified that values for A of up to 210 may be 
excluded, if mSUGRA physics remains unobserved at the LHC with ^/s = 14 TeV 
and a luminosity of 10/6~^. The study included the dominant one loop corrections 
to the Higgs potential. Having identified the magnitude of errors associated with one 
loop fine tuning calculations in this work, and noting that fine tuning may be slightly 
reduced by varying A^ and tan/3, the limits on A found here are consistent with this 
previous work. 

For the LHC to operate at ^/s = 10 TeV and obtain a luminosity of 100 pb^^, 
gluinos of up to 600 GeV may be discovered. However, if gluino masses are excluded 
to this level, A is only excluded to a value of 14. As has been discussed, charginos 
and neutralinos can be quite light, but their signal events are difficult for LHC to 
extract from the background, owing in part to a decreasing M^ — M^ mass gap 
as |yu| decreases [321 SO]- An Atlas study [H] of the trilepton signal from chargino- 
neutralino production found that 30/6^^ luminosity at 14 TeV is needed for a Sex 
discovery significance for M2 < 300 GeV and /i < 250 GeV [12] . 

2.4 Summary 

Supersymmetry was introduced to solve the electroweak hierarchy problem and to 
avoid the large fine-tuning in the SM Higgs sector associated with the Planck or 
gauge unification scale when quantum corrections are included. While this hierarchy 
problem is solved by TeV-scale supersymmetry, the non-observation, so far, of SUSY 
states means that the MSSM has acquired some residual amount of fine-tuning related 
to unnatural cancellations in the SUSY breaking sector. 

The fine tuning measure A provides a quantitative test of SUSY as a solution to the 
hierarchy problem and measures the "tension" required to satisfy the scalar potential 
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minimum condition t>^ ~ —ni^^^y/X, for a combination of soft masses m^u^y ~ TeV, 
with an effective quartic coupling A remaining perturbative and v ~ (9(100) GeV. 
Although the exact upper limit on the fine tuning beyond which a theory fails to 
solve the hierarchy problem is debatable, it is preferable for a given model to have 
a parameter space configuration corresponding to the lowest value of A. The fine 
tuning measure has been evaluated here at two-loop order with particular attention 
given to threshold corrections and the tan /3 radiative dependence on the parameters. 
Such effects on fine tuning were not fully considered in the past and turned out to 
reduce fine tuning significantly. 

The determination of the fine tuning measure for mSUGRA included the the- 
oretical constraints (radiative EWSB, avoiding charge and colour breaking vacua), 
and also the experimental constraints (bounds on superpartner masses, electroweak 
precision data, 6— )'S7, 5^— T-yUyU and muon anomalous magnetic moment). 

The analysis found that a fine tuning of 1 part in 0(10) remains consistent with the 
current experimental constraints. This is much more promising than the commonly 
reported limit of 1 part in O(IOO) frequently repeated in the literature based on very 
primitive calculations. Moreover, these results indicate a preference for a Higgs mass 
at ~ 115 GeV, which was found without application of the LEPII Higgs mass bound. 
This may be regarded as a theoretical prediction for the Higgs mass based on the 
principle of naturalness. Chapter 4 of this thesis will also discuss how non-minimal 
theories may improve naturalness beyond the degree found possible in mSUGRA. 

The spectrum corresponding to the minimum value of the fine tuning shows sim- 
ilarities to the SPS2 scenario with light neutralinos, charginos and gluinos (corre- 
sponding to light /i, mi/2) and heavy squarks and sleptons corresponding to large mo, 
near the focus point limiting value |l3l HI]. It provides the "best" estimate for the 
SUSY spectrum given the present experimental bounds. 

Increasing ruh above the minimum fine tuned value causes A to increase expo- 
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nentially fast. One obtains A = 100 (1000) for a scalar mass rrih = 121 (126) GeV, 
respectively. Ultimately the question whether the SUSY solution to the hierarchy 
problem has been experimentally tested relies on what value of fine tuning represents 
the limit of acceptability. Given a value, one can determine the range of parameter 
space that is still acceptable. For the case that the fine tuning measure should sat- 
isfy A < 100, the corresponding ranges for the superpartners masses and mSUGRA 
parameters values have been determined to be relevant for SUSY searches. 



Chapter 3 



Fine tuning of dark matter within 
minimal supergravity 



In order to test a theory for naturalness, it is not sufficient to test whether just one 
observable, such as the electroweak scale, is fine tuned. One must test all observables. 
If any observable is found to require a large degree of fine tuning, then the theory must 
be deemed to be "un-natural." Due to the electroweak scale being notably vulnerable 
to fine tuning issues, it is sensible to first test for naturalness of electroweak symmetry 
breaking (EWSB). Having found in the previous chapter fine tunings of less than 1 
part in 10 consistent with current observations, attention is now turned to testing 
other physics within the mSUGRA model. 

An exhaustive study of naturalness for any given model would investigate many 
different sectors of physics, such as fiavour physics, QCD phenomenology and cos- 
mology (infiation). For this thesis, the scope is restricted for consideration of the 
electroweak sector. As mentioned before, this is of particular interest for possible 
imminent discoveries. The phenomenology of other sectors are also not necessarily 
connected with electroweak physics, allowing this separation of analyses. However, 
as alluded to in earlier chapters, there are appealing motivations for unifying physics. 
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and so unified models would need to check for naturalness in observables of all con- 
nected sectors. 

In addition to motivating supersymmetry as a solution to the hierarchy problem, 
the presence of a dark matter candidate was argued. Evidence for dark matter is 
obtained from analysis of galactic rotation curves. The expected rotation from the 
observed luminous matter using standard gravity does not match the results found. 
The existence of diffuse halos of matter that do not interact significantly with light, 
ie "dark" matter, are then required to match the observations. 

Furthermore, the power spectrum of the cosmic microwave background exhibits 
peaks from which the ratio of luminous matter to total matter content of the universe 
can be inferred. The WMAP experiment [17] determines this ratio to be ~ 0.2. 
The majority of matter is then non-luminous. Simulations of galaxy formation also 
suggest that the dominant contribution to dark matter is non-relativistic, or "cold" 
dark matter. 

The choice of R-parity within the MSSM stabilised the lightest supersymmetric 
particle (LSP) with respect to decay. The spectrum of mSUGRA is presented in 



Section 13.11 leading to identification of the LSP as having electroweak interactions 
(the gravitino is assumed heavy). Thus, as well as testing for naturalness of EWSB 
on introducing SUSY, one must also test the physics of this dark matter candidate 
to see if the electroweak sector is natural. 

If R-parity is not realised, making the LSP unstable to decay, or the LSP is actually 
not related to the electroweak sector (such as an axino, gravitino, or a hidden sector 
state), then this analysis is not relevant. However, if the dark matter candidate 
considered here provides any thermal contribution to the relic energy density, the 
restriction on parameter space and naturalness constraints presented are applicable. 



Section 13.21 discusses how the thermal relic energy density of a stable particle can 
be calculated. Section 13.31 reviews the constraints obtained on dark matter physics 
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from direct detection searches. Section 13.41 presents the allowed parameter space 



found on imposing the dark matter constraints, and Section [331 then presents the fine 
tuning of the relic density in the experimentally allowed parameter space. 



3.1 MSSM spectrum 



The spectrum of the superpartner states is controlled by the unknown effective soft 
SUSY breaking terms, and the renormalisation group equations (RGEs). For the case 
of mSUGRA, universal soft SUSY breaking parameters are fixed at the gauge unifi- 
cation scale. The MSSM renormalisation group equations |1H] , which are dominated 
by the third generation couplings, then drive the stops and staus to be the lightest 
squarks and sleptons. 

As the 1-loop RGEs (see Appendix 1X1) obey the relation, dlnM^/dt = dlng'^/dt, 
the gaugino soft mass, M^, run approximately equal to the gauge couphngs squared, 
g"^. The universality of gaugino masses is broken on running away from the UV scale, 
but the ratio Mi : M2 : M3 then roughly follows a^ : ag • '^3- ^^ ^^^ electroweak 
scale, this ratio is approximately given by 1 : 2 : 6 (using an SU{5) normalisation 
for hypercharge) , and in terms of the mSUGRA soft mass parameter, fn^,^, which 
sets the gaugino mass at the gauge unification scale, Mi{m^) ~ OAm^,^- The ratio 
relations are modified by threshold corrections and higher loop effects. 

The neutralino mass matrix in the {B^, W^, h\, h^ basis is given by. 
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(3.1) 



where Sa, (c^) = sin/3 (cos/3) and s^y, {c^/) = sm9^/ (cos6'p,/). B^ (W) is the super- 



partner of the hypercharge (weak) gauge boson and known as the Bino (Wino), and 
hi 2 are the Higgsinos. Note that the gluino cannot mix with these states as colour 
is an unbroken non-Abehan symmetry. The chargino mass matrix in the {W^,h^) 
basis is, 
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(3.2) 



The magnitude of fi is fixed to reproduce the correct electroweak gauge boson 
masses. The vacuum minimisation equations of the Higgs potential tend to lead 
to the hierachy, n > M2 > Mi. However, in the region of parameter space where 
the scalar focus point exists, one finds M2 > /U > Mi. In both cases, the lightest 
neutralino is lighter than the charginos. Neither the magnitude of fi nor 171^,2 can be 
too small, as experimental limits on the superpartner masses forbid such cases. For 
mSUGRA models, these limits disallow fi ^ M^. 

For the common hierarchy, n > M2 > Mi, the lightest neutralino, Xi5 is domi- 
nantly Bino-like and X2iXza ^^^ Wino-like and Higgsino-like respectively. Similarly, 
the light (heavy) chargino is Wino-like (Higgsino-like). The ratio of the two lightest 
neutralinos and gluino then also roughly follows the ratio, «! : ^2 • '^3- 

For the scalar focus point region, M2 > /U ^ Mi, there is significant mixing of the 
gaugino/Higgsino fields in the neutralino and chargino mass eigenstates. In the limit 
Ml ~ yU, and large tan/3, the lightest neutralino mass is approximately given by. 
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(3.3) 



The states {xi 1 X2 3} ^-re all 0{Mi), nearly degenerate with the lightest neutralino, 
and {X2'!X4} ^'^^ Wino-like with mass 0(M2 ~ 2M^). For the smallest experimen- 
tally allowed values of M^ in mSUGRA, the ratio of the gluino mass to the lightest 



neutralino increases from ~ 6 to 0(10). As the coloured gluino is one of the most 
promising particles for early detection at the LHC (assuming a light mass), this mass 
hierarchy will be referred to in the following results which plot the lightest neutralino 
mass. The mass ratios within the spectrum will be important to understand which 
regions of parameter space are detectable in the first runs of the LHC. 

In order for the LSP to be a dark matter candidate, the regions of parameter 
space where the LSP is electrically charged (eg stau or stop) must be ruled out. The 
stop mass squared matrix is given by. 



m? + mf + L m^ [A^ — fi cot f3) 
> m^ {A^ — ji cot (3) m? + ml + R^ 



(3.4) 



where L^ = (| — | sin^ 6^r^ m| cos 2/3, and -R^ = | sin^ ^w^z cos 2/3. 

There is a similar structure for sbottoms and staus (although with different tan/3 
dependence and distinct L^, R^ definitions). It is usually the mass mixing induced by 
large Aq/itlo that leads to a small stop mass eigenstate. In the analysis of fine tuning 
for the electroweak scale, small A^ magnitudes were preferred on naturalness grounds, 
and so large sfermion mixing tends to require large fine tuning. As we are interested in 
determining the low fine tuned regions, this constraint is not problematic. Assuming 
that the gravitino is heavy, the lightest neutralino then is identified as the LSP. This 
particle will interact through the electroweak interactions. 

3.2 Thermal relic energy density 

In the Big-Bang scenario, an abundance of matter in the present day may result from 
an asymmetry between matter and antimatter. Alternatively, stable particles may 
leave thermal equilibrium with its annihilation products, "freeze-out", to generate a 
relic density in the present day. The latter case is now discussed for its predictions 
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of the relic density. The freeze-out will occur roughly when the thermally averaged 
annihilation rate of this particle becomes slower than the expansion rate of the uni- 
verse: 

n{a,,,,,v) r^ H (3.5) 

where n is the number density of the particle freezing out, (cr^nnV) is the thermally 
averaged annihilation cross section weighted by the relative velocity, and H is the 
Hubble constant. 

The ensemble of frozen-out particles then expand isenthropically, allowing the 
present number of particles in a comoving volume to be calculated. The current 
contribution to the energy density of the universe can then also be determined. If the 
particles are non-relativistic at the time of freeze-out, the thermal average is given 
by: 

POO 

{(T^nnv) ^ f^ / dvv'ia^^^^v)e-^^'/' (3.6) 

where x = ni-^j^/T and m-^-^ is the mass of the stable particle (dark matter). The 
annihilation cross section can be calculated perturbatively, however for slow moving 
dark matter there can be important non-perturbative effects due to interactions before 
annihilation. This is known as the Sommerfeld effect. For further discussion of 
this effect, see ^9] and the contained references. It turns out that for electroweak 
interactions the Sommerfeld effect is only important for precision calculations at a 
level beyond the error introduced from the following approximations. 

The temperature at which the freeze-out occurs is determined from eq (13. 5p . and 
is parametrised by the value of x at freeze-out |50] : 

_ J^ / [3.85X1017 GcV] gDM'^DM('^^nn^).^ \ , . 
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where g^^^ is the number of degrees of freedom of the dark matter candidate. The 
g^^s parameter is the effective number of relativistic degrees of freedom contributing 
to the entropy density of the universe. For T ~ 300 GeV, (7*^ = 106.75 (assuming 
only Standard Model particles contributing), and at T = 1 GeV, g^:s ~ 80. The 
contribution to the energy density is then given by: 



n^^e = 'T \T~ ^here J{xf) = / dxx~-'{a^^^v) (3. 



where the approximation that all annihilations cease after "freeze-out" has not been 
applied. The "annihilation integral," J{xp), accounts for the reduction in particle 
number after "freeze-out." The thermal relic will contribute to the non-baryonic mat- 
ter energy density observed in the universe today. The 5- year data from WMAP [17] 
suggests that ^n.b.m.^'^ = 0.1099 ± 0.0062, with the la error stated. 

The thermal relic abundance increases as the annihilation cross section is reduced. 
The unitarity limit places an upper bound on dark matter masses with such thermal 
relic density at O(IO^) TeV [51]. For dark matter with electroweak interactions, a 
mass of O(TeV) is required to obtain the observed energy density. Supersymmetry 
provides such a dark matter candidate as the lightest superpartner state, being sta- 
bilised by an R-parity. Note that the mass scale is the same as that argued for solving 
the hierarchy problem. 

3.3 Direct dark matter searches 

In order to verify the existence of dark matter without inference from galactic rota- 
tion curves and cosmological studies, there are several types of experiment ongoing 
for this purpose. The experiments are either concerned with direct or indirect detec- 
tion of dark matter. Indirect detection aims to observe the annihilation products of 
dark matter from astrophysical sources. The constraints from such experiments are 
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currently subdominant to direct detection experiments for mSUGRA dark matter. 
Interpretation of the results also relies on astropliysical models and precise catalogs 
of background sources leading to large uncertainties. These types of experiment will 
not be discussed further. 

Direct detection of dark matter may be acheived by measuring recoils of nuclear 
matter from collisions with dark matter passing through the experiment. From the 
rotation curve of our galaxy, the expected local dark matter density is 0.3 GeV cm~^ 
with an 0(2) factor uncertainty (see |53j and contained references). The large uncer- 
tainty is a consequence of the difficulty determining the distance of our solar system 
from the galactic centre (our angular velocity is well known though). 

Computational simulations suggest that the structure of observed galaxies is con- 
sistent with a smooth distribution of dark matter on the scales that are sampled by 
the direct detection experiments. The relative velocity, f , of the Earth to the dark 
matter halo is believed to be 220 kni/s varying by 0(10) km/s seasonally, and as 
direct dark matter detection experiments run over a timescale of years, this covers a 
O(10"^Parsec) distance. 

If the calculated mSUGRA dark matter candidate does not saturate the non- 
baryonic matter energy density as determined by WMAP, the missing energy density 
may be accounted for by an alternative dark matter state. The identity of alternative 
dark matter will not be discussed here, however such a situation affects how the 
constraints from direct detection dark matter experiments must be handled. 

In the following analysis, the non-baryonic matter energy density is always as- 
sumed to match the result determined by WMAP. The local density of the mSUGRA 
dark matter, p , is then assumed to scale proportionally as the fraction of mSUGRA 
dark matter to total dark matter present in the universe does. This is reasonable for 
(approximately) coUisionless dark matter, which numerical simulations suggest is a 
property of dark matter. 
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The differential scattering rate per unit recoil energy and unit target mass is 

/ aF\E^)£l^d\ (3.9) 



dR ^ Px 



dE ^^^ 



where R is the number of recoil events measured, Eji is the nuclear recoil energy, a is 
the elastic LSP-nucleon interaction cross section (for which an exclusion limit is to be 
found), /i^N is the reduced mass of LSP-nucleon system, m^ is the LSP mass, F{q) is 
the nuclear form factor, v^i^ is the minimum velocity required to give nuclear recoil 
energy Eji, v^^c is the escape velocity of dark matter in our galaxy (boosted into Earth 
frame), and f{v) is the velocity distribution of LSP in the local neighbourhood. 

The non-observation of recoil events above the background rate places an upper 
limit on the LSP-nucleon cross section. The elastic scattering amplitudes can be 
separated into the contributions that are proportional to the spins of the scattering 
states (spin dependent), and those that are spin independent. As the two contribu- 
tions do not interfere, separate limits can be placed on the spin dependent and spin 
independent cross sections. 

It is the spin independent cross section limit that is the dominant constraint for 
neutralino dark matter. This may be expected as the spin independent amplitudes 
add coherently, whereas the spin dependent amplitudes are incoherent. In order to 
have large nuclear recoils, heavy nuclei are preferred for the experimental set up, and 
so the large number of nucleons lead to a relatively larger cross section for the spin 
independent case. 

Figs 13.11 and 13.21 show the experimental limits produced by the CDMS II experi- 
ment [52]. This is currently the strongest constraint on the LSP-nucleon scattering 
cross section for elastic dark matter. The points plotted correspond to points in the 
mSUGRA parameter space that satisfy the experimental limits listed in Table 12.11 as 
well as the relic density constraint. 
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Figure 3.1: Spin independent cross section for LSP-proton scattering, with an unrestricted 
Higgs mass. The sohd line is the CDMS II hmit, which is only applied in (b) and (d). 
The points satisfy the constraints in Table [2TT] and have Qh'^ < 0.1285. Points with darker 
shading are within 3a of the total non-baryonic matter relic density, O/i^ = 0.1099 it 3 x 
0.0062. The labels distinguish between the dominant LSP annihilation mechanisms. The 
h^ {H^,A^) resonance is significant for region 1 (5) with points coloured red (black). Region 
3 (4) realises stau (stop) co-annihilation with points coloured green (blue), and region 2, 
coloured purple, has increased Higgsino components. 

The figures show that the direct dark matter constraint does exclude certain parts 
of the parameter space, whether the Higgs mass limit is applied or not. The phe- 
nomenology of the labelled regions is discussed in the following section. As future 
direct dark matter detection experiments will improve the reach by orders of mag- 
nitude in the next few years, there is competition with the collider experiments for 
discovery of supersymmetric states. 

The procedure of scaling the local LSP density by the relic density fraction means 
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Figure 3.2: Spin independent cross section for LSP-proton scattering, with ruh > 
114.4 GeV. The solid line is the CDMS II limit, which is only applied in (b) and (d). The 
points satisfy the constraints in Table [2T] and have Qh"^ < 0.1285. Points with darker shad- 
ing are within 3a of the total non-baryonic matter relic density, Qh'^ = 0.1099 it 3 x 0.0062 
(for all following plots also). The labels distinguish between the dominant LSP annihilation 
mechanisms. The h^ {H^,A^) resonance is significant for region 1 (5) with points coloured 
red (black). Region 3 (4) realises stau (stop) co-annihilation with points coloured green 
(blue), and region 2, coloured purple, has increased Higgsino components. 

that a weaker limit is obtained for such points. This is the reason for finding some 
of the lighter shaded points (those that do not saturate the relic density) above the 
CDMS II line in subfigures (b) and (d), in which the constraint is appled - the line 
assumes that the LSP saturates the local dark matter. 

The dominant astrophysical uncertainty is generated by the uncertainty in the 
local dark matter density, which is neglected for the limits produced by experimental 
groups. It is standard practice to assume the best fit values for the astrophysical 
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variables, and so the experimental limits shown in the plots should be considered to 
vary by an 0(2) factor according to the astrophysical input. 

3.4 Phenomenology of allowed parameter space 

In addition to the constraints on the parameter space applied in the previous chapter 
(as listed in Table l2TT]) . the thermal relic density of the neutralino dark matter is now 
restricted such that, Qh"^ < 0.1285, for agreement with the WMAP 5 year results but 
allowing a 3o" error deviation. The CDMS II limit is also applied as detailed in the 
previous section. The dark matter relic abundance and LSP-nucleon cross sections 
are calculated using MicrOMEGAs 2.2 [54j. 

The allowed parameter space of mSUGRA is strongly restricted by the condition 
of requiring a small dark matter relic density. This is a consequence of the neutralino 
LSP being dominantly Bino-like for most of the parameter space allowed before apply- 
ing the dark matter constraint, as mentioned in Section [XT! The small hypercharge 
coupling leads to a small annihilation cross section and, as Section [3]2] demonstrated, 
the thermal relic density is inversely proportional to the annihilation cross section. 

The mechanisms for generating an increased annihilation cross section, and so 
a reduced relic density, distinguishes the allowed parts of parameter space that re- 
main after applying the dark matter constraints. The universal structure of the soft 
SUSY breaking terms also leads to specific predictions of the mSUGRA spectrum 
that controls this phenomenology. In a more general supersymmetric model, the 
phenomenology is less restricted on application of a fixed number of constraints. 

Mechanisms that can reduce the relic abundance include the presence of co- 
annihilation channels, where states nearly degenerate with the LSP interact with 
the LSP to produce final states without any superpartner states. If the parameters 
that lead to nearly degenerate states are distinct, this is expected to require fine tun- 
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ing. Another inechanism is mixing with Higgsino states which is discussed further in 
Section [3A21 

The final mechanism to be discussed involves resonant annihilation. As the neu- 
tralino dark matter is non-relativistic at the time of freeze-out, the centre of mass 
energy for annihilation processes roughly obeys the relation, ^/s ~ 2m o- If the 2m o 
happens to be equal a bosonic state (such as a Higgs field), then the relic abundance 
will be depleted by the resonance. Again, this scenario suggests some degree of fine 
tuning is required. 

The following subsections discuss the specific phenomenology of certain mech- 
anisms that are found to lead to the relic density constraint being satisfied. The 
remaining allowed parameter space is separated into five regions, and the numbering 
is used to refer to each region in the plots of this chapter. 

The case of a (left-handed) sneutrino LSP has been overlooked so far in the discus- 
sion. This is a viable dark matter candidate, but the large electroweak interactions 
lead to direct detection experiments excluding the candidate. The limits on super- 
partner masses and UV structure of mSUGRA also restricts the region of parameter 
space where the sneutrino is the LSP to a relatively small region, when the direct 
search limits are ignored. This situation will not be discussed further. 

In past literature, the "bulk" region has also been identified as consistent with the 
relic abundance. In this region, the Bino annihilation is increased by the exchange 
of light sleptons in the t and u-channel, which occurs at small rn-^,^ and small m^. 
Figs 13.51 and 13.61 demonstrate that this region is now excluded, which is a consequence 
of the improved experimental constraints on superpartner masses. 

3.4.1 Region 1: h^ resonant annihilation 

The lower bound for lightest neutralino mass in the mSUGRA scenario is 46 GeV 
for most of the parameter space. The centre of mass energy for LSP annihilation 
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Figure 3.3: Plot of LSP mass vs LSP composition in gaugino-Higgsino basis, with the 
Higgs mass unrestricted. The overlap with respect to the Bino is plotted in (a), the Wino 
in (b), the down- type Higgsino in (c) and the up- type Higgsino in (d). The labels dis- 
tinguish between the dominant LSP annihilation mechanisms. The h^ {H^,A^) resonance 
is significant for region 1 (5) with points coloured red (black). Region 3 (4) realises stau 
(stop) co-annihilation with points coloured green (blue), and region 2, coloured purple, has 
increased Higgsino components. 

then obeys, i/i > 92 GeV. For all but the smallest values of tan/3, the light Higgs 
resonance will be accessible in some region of the parameter space allowed before 
dark matter constraints. The width of the h^ resonance is 0(3 GeV), and so it may 
not require significant amounts of fine tuning. As for all the regions identified in this 
section, the degree of fine tuning will need to be tested. 

Figs 13.31 and 13.41 demonstrate that only a small range of neutralino mass (50 
to 62 GeV) actually leads to a satisfactory relic abundance, and the LSP is indeed 
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Figure 3.4: Plot of LSP mass vs LSP composition in gaugino-Higgsino basis, with mh > 
114.4 GeV. The overlap with respect to the Bino is plotted in (a), the Wino in (b), the down- 
type Higgsino in (c) and the up- type Higgsino in (d). The labels distinguish between the 
dominant LSP annihilation mechanisms. The h^ (H^,A^) resonance is significant for region 
1 (5) with points coloured red (black). Region 3 (4) realises stau (stop) co-annihilation 
with points coloured green (blue), and region 2, coloured purple, has increased Higgsino 
components. 



Bino-like. The Higgs mass constraint also does not considerably restrict this region of 
allowed parameter space. A light Higgs resonance at this energy scale will dominantly 
decay to bb products. It is then larger tan /3 that will lead to the large bottom Yukawa 
couplings to make the resonance most efficient. As the light Higgs is CP-even, the 
annihilation will dominantly proceed through p-wave annihilation. This is then less 
effective as other methods for reducing the relic density. 



Figs 13.11 and 13.21 show that the spin independent cross section for this region of 
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parameter space is significantly below the current experimental limits so the most 
promising hope for early discovery is at colliders. The gluino is roughly six times 
heavier than the LSP in this region, and as the LSP (and gluino) here is only just 
above the current experimental limit, it will not require reaching the design centre of 
mass energies at colliders to generate these states in significant numbers. 

Figs 13.51 and 13.61 indicate that this region is around the scalar focus point region, 
with a large m^jm^,^ ratio. The region also extends to large A^ parameters. This free 
variable can be used to introduce appropriate amounts of stop mixing to modify the 
Higgs quartic couplings and so shift the Higgs mass to be on resonance. As observed 
in the last chapter, this requires a sacrifice in fine tuning of the electroweak scale and 
so such limits for this region are not desirable on naturalness grounds. 

3.4.2 Region 2: Bino-Higgsino mixing 

On the edge of parameter space where electroweak symmetry breaking is allowed, the 
SUSY yU parameter is driven towards zero by the Higgs potential vacuum minimisation 



conditions. Section 13.11 indicated that the superpartner mass limits restrict the [i 
parameter from being much smaller than M^. In this region, there is a significant 
Higgsino component within the LSP. 

The Higgsino can annihilate via a Higgsino t-channel to produce electroweak gauge 
bosons. The couplings are relatively large and so the annihilation is efficient. In fact, 
for pure Higgsino states, the relic abundance would be generally far below the WMAP 
result. An intermediate mixture of Bino and Higgsino is preferred for generating the 
observed relic abundance. 



Section 13.11 also mentioned that there is a near degeneracy in this region with 
the light chargino and all but the heaviest neutralino. There will therefore also be 
coannihilation channels from the x? Xi ^^ ^-^^d Xi X2 ^ vertices. 

Figs 13.31 and 13.41 demonstrate that in this region, the LSP is 60 to 75% Bino, 10 
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Figure 3.5: Plots of the allowed mSUGRA parameter space subject to all constraints 
listed in Table [2m and the dark matter constraints, and in (b) spin-dependent cross section 
vs LSP mass (where the experimental limit is off the top of the plot). The Higgs mass 
is unrestricted. The stepped choices of tan/3 is responsible for the lines of points in (e) 
and (f). The labels distinguish between the dominant LSP annihilation mechanisms. The 
h^ {H^, A^) resonance is significant for region 1 (5) with points coloured red (black). Region 
3 (4) realises stau (stop) co-annihilation with points coloured green (blue), and region 2, 
coloured purple, has increased Higgsino components. 
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Figure 3.6: Plots of the allowed mSUGRA parameter space subject to all constraints 
listed in Table [2m and the dark matter constraints, and in (b) spin-dependent cross section 
vs LSP mass (where the experimental limit is off the top of the plot). The constraint, 
nih > 114.4 GeV, is also applied. The stepped choices of tan/3 is responsible for the lines 
of points in (e) and (f). The labels distinguish between the dominant LSP annihilation 
mechanisms. The h^ {H^,A^) resonance is significant for region 1 (5) with points coloured 
red (black). Region 3 (4) realises stau (stop) co-annihilation with points coloured green 
(blue), and region 2, coloured purple, has increased Higgsino components. 
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to ~20% of each Higgsino type, and up to 5% Wino when the non-baryonic matter 
rehc density is saturated by the LSP contribution. In the case that the rehc density 
is sub- saturating, the LSP may be half each Higgsino type (and < 1% Bino). 

The Higgs mass constraint does not affect this region at all, but Figs 13.11 and 



13.21 show that the CDMS II constraint is restrictive for light neutralino masses. The 
large Higgsino components lead to large couplings with nucleons. It is this region of 
parameter space that the direct detection experiments are currently probing. 

3.4.3 Region 3: Stau co-annihilation 

The stau is driven relatively light in the region of parameter space where m^i^/mQ is 
large as observed in Figs 13.51 and 13.61 This is expected, as the diagonal entries in the 
stau mass squared matrix are 0{m1). 

Figs 13.31 and 13.41 demonstrate that the LSP mass is correlated with the Wino and 
Higgsino components, being lighter for increased non-Bino components. This is just 
a consequence of the mSUGRA pattern for the neutralino mass matrix affecting the 
mixing as the magnitude of m^ /2 is varied. The Higgs mass constraint signifcantly 
restricts this region of parameter space, as it is mostly found for low tan/3 values 
where the Higgs mass tends is light. 



Figs 13.11 and 13.21 show some restriction is provided by the CDMS II constraint 
when the Higgs mass is ignored. However, with the Higgs mass constraint, the points 
in this region are over an order of magnitude below the direct detection limits. The 
variation in the spin independent cross section for varying neutralino mass does have 
some kinematic origin, but the dominant influence is the varying Higgsino component 
which controls the effectiveness of coupling to nucleons. 
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3.4.4 Region 4: Stop co-annihilation 

The stop co-annihilation region is in a similar region as the stau co-annihilation region. 
However, as indicated in Figs l3.5l and l3.6[ a large Aq parameter is necessary to induce 
large stop mixing. The mixing was less effective for staus due to the smaller tau 
Yukawa coupling. 

Figs 13.31 and 13.41 demonstrate that the non-Bino components of the LSP are ex- 
tremely small (< 0.5%), and so the spin independent cross sections plotted in Figs 13. II 
and l3.2l are therefore relatively small. These points are also found for low to moderate 
tan/3, and so the Higgs mass constraint significantly restricts this region. At large 
tan P the off-diagonal terms in the stop mass matrix are reduced, and so there is less 
mixing. 

3.4.5 Region 5: H'^,A^ resonant annihilation 

Similarly to region 1, where there is resonance with the light Higgs state, there exists 
a region of allowed parameter space where resonance with the heavy Higgs states 
aids in reducing the relic density. Due to the near degeneracy of these states and 
the 0(30 GeV) widths, it is not possible for this parameter space scan to resolve each 
resonance separately. 

Figs 13.31 and 13.41 demonstrate that this region has a similar LSP composition to 
the stau co-annihilation region. The dependence of the spin independent cross section 



shown in Figs 13.11 and 13.21 then also has a similar behaviour. 

However, the structure of the mSUGRA mass spectrum requires large m^ ^ and 
large tan /3 in order to find this resonance as indicated in Figs 13.51 and 13.61 The 
annihilation via the CP odd A^ can proceed via s-wave annihilation, but as for the 
h^ resonance, the CP-even H^ resonance must proceed via p-wave annihilation. A 
larger suppression is generally found when closer to the A^ resonance. 
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3.5 Fine tuning results from parameter space scan 

To quantify the fine tuning of the rehc density, A^, the Barbieri-Giudice measure is 
again used, 



A = max 



dlnUh^ 



dlnq 



(3.10) 



except that the parameters tested now include q = {wiq, m^ ,2, Aq, tan /?}. The Ciafaloni- 
Strumia measure was considered for the measured parameters, {171^,771^}, but always 
found to be subdominant and so effectively ignored. The set of parameters are not 
all considered to be fundamental though. This compromise was necessary for using 
the numerical MicrOMEGAs package in order to determine the relic density. A routine 
for evaluation of fine tuning is also not built in the package as the fine tuning for the 
electroweak scale is within SOFTSUSY. 

The procedure to evaluate the relic density fine tuning then used a finite difference 
method, checking the variation under small finite changes in the input parameters. 
Convergence tests were performed to confirm that the results were reliable to an 
0(1%) level. The theoretical uncertainty for most parts of the parameter space is 
expected to be at this level given comparisons with DarkSUSY. For earlier analyses of 
the fine tuning of the relic density, see [SS, \5E\ and the references therein. A similar 
approach to that adopted here was used in those studies, although their scope did not 
cover a full parameter space scan, nor use the now available experimental constraints. 

The following results demonstrate the degree of fine tunings present on application 
of the dark matter constraints. The results for fine tuning of the relic density and the 
electroweak scale are both presented, as well as plots to demonstrate the worst fine 
tuning present for a given point in the parameter space. For the latter, the worst fine 
tuning is defined to be max[A, A^], following a similar procedure as A being defined 
as max[{Ap}]. 
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This procedure of combining fine tunings for different observables is not rigourously 
motivated. It was argued in the introductory chapter that the fine tuning of any one 
observable is not well defined, with many different approaches equally reasonable. 
One ambiguity was in choosing an appropriate parameterisation of the fundamental 
parameters for a given observable (dependent on its mass dimension). The "worst" 
fine tuning is considered here for a straightforward interpretation of the results - the 
measure indicates the greatest degree of fine tuning necessary for either considered 
observable to assume its physical value. 

3.5.1 Dependence on dark matter constraints 

In order to compare with the results of the previous chapter, Fig 13.71 reproduces its 
results with the relic density constraint and subsequently the CDMS II constraint 
applied. The solid line and shaded area indicate the filling obtained without the dark 
matter constraint. The points (for which the dark matter constraint is applied) fill 
up only a fraction of that space. 

Region 1 {h^ resonance) is observed to occupy a region of low fine tuning around 
the LEP II Higgs mass bound. For points that sub-saturate the non-baryonic relic 
density, the smallest fine tuning is accessible. In fact, with the CDMS II constraint, 
the global minimum for the parameter space scan in the fine tuning of the electroweak 
scale and the overall fine tuning is found at the LEP II limit (just as was the case 
when ignoring the dark matter constraints) . The impact of including the fine tuning 
of the relic density leads to slightly increased fine tuning in this limit and so both 
the electroweak scale and relic density are fine tuned to roughly 1 part in 15 at this 
minimum. 

Region 2 (Bino-Higgsino mixing) is found to fill in the large Higgs mass region. 
The fine tuning of the relic density is greater than that for the electroweak scale at the 
lowest fine tuned points (low Higgs mass). However, these points are also excluded 
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Figure 3.7: Plots of fine tunings vs Higgs mass. The solid line and shading indicates the 
region allowed without any dark matter constraints. No direct detection experimental limits 
are applied in (a) and (c), but in (b) and (d) the CDMS II limit is applied for the points 
plotted. The fine tuning of the electroweak scale is plotted in (a) and (b), and the maximum 
fine tuning of the electroweak scale and relic density is plotted in (c) and (d). The labels 
distinguish between the dominant LSP annihilation mechanisms. The h^ {H^, A^) resonance 
is significant for region 1 (5) with points coloured red (black). Region 3 (4) realises stau 
(stop) co-annihilation with points coloured green (blue), and region 2, coloured purple, has 
increased Higgsino components. 



by the CDMS II constraint. For r/i^ > 118 GeV, the influence of dark constraints on 
the line of mininum fine tuning is negligible and the relic density fine tuning is not 
dominant in this limit except for the largest Higgs masses. 

Region 3 (Stau coannihilation) seems to have strips of allowed space at low Higgs 
mass. This is in fact a result of the discrete stepping in tan /3 for the parameter space 
scan. A smoother scan across tan/3 fills in the gaps between the 'strips.' It is observed 
that region 2 covers a wide range of Higgs mass. The relic density fine tuning does 
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not noticeably affect the overall fine tuning. The impact of the CDMS II constraint 
is to eliminate points with low Higgs mass, at low tan/3. These points are far below 
the LEP II bound anyway. 

Region 4 (Stop coannihilation) possesses large fine tuning of the electroweak scale, 
A > 200, which is frequently increased further by the relic density fine density at low 
tan p. This was expected from the large A^ parameters required to generate the stop 
mixing. 

Region 5 {H^, A^ resonance) also requires moderately large fine tuning, with the 
relic density fine tuning regularly dominating. The light Higgs mass ranges from just 
below the LEP II bound to 120 GeV. 

3.5.2 Dependence on SUSY Lagrangian parameters 

The following results now always apply both the relic density and CDMS II constraints 
for the points plotted. The relic density fine tuning is also independently plotted for 
comparison and analysis. The plots of fine tuning for the electroweak scale and overall 
fine tuning indicate the region allowed when ignoring the dark matter constraints. 

Fig 13.81 shows the dependence of fine tunings on the scalar soft mass. The impact 
of the dark matter constraints is to remove a large part of the parameter space with 
intermediate m-^^,2 with smallest fine tuning (when the Higgs mass is unrestricted). 

Region 1 {h^ resonance) is seen to populate the 'dip' of fine tuning in rriQ found on 
application of the LEP II bound. The relic density fine tuning also extends to very 
small values. Region 2 (Bino-Higgsino mixing) is seen to have a well-defined structure 
in the fine tuning of the electroweak scale, but the relic density fine tuning frequently 
raises the overall fine tuning significantly. Region 3 (Stau coannihilation) includes 
some very low relic density fine tuned points, but these are eliminated by the Higgs 
mass constraint. Region 4 (Stop coannihilation) and Region 5 {H^, A^ resonance) fill 
in the space between the other regions. 
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Figure 3.8: Plots of fine tunings vs UV scalar soft mass, tjiq. The fine tuning of the elec- 
troweak scale is plotted in (a) and (b), the fine tuning of the relic density in (e) and (f), and 
the maximum fine tuning between the two in (c) and (d). The labels distinguish between 
the dominant LSP annihilation mechanisms. The h^ (H^,A^) resonance is significant for re- 
gion 1 (5) with points coloured red (black). Region 3 (4) realises stau (stop) co-annihilation 
with points coloured green (blue), and region 2, coloured purple, has increased Higgsino 
components. 
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Figure 3.9: Plots of fine tunings vs UV gaugino soft mass, mi2. The fine tuning of the elec- 
troweak scale is plotted in (a) and (b), the fine tuning of the relic density in (e) and (f), and 
the maximum fine tuning between the two in (c) and (d). The labels distinguish between 
the dominant LSP annihilation mechanisms. The h^ {H^,A^) resonance is significant for re- 
gion 1 (5) with points coloured red (black). Region 3 (4) realises stau (stop) co-annihilation 
with points coloured green (blue), and region 2, coloured purple, has increased Higgsino 
components. 
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Figure 3.10: Plots of fine tunings vs UV trilinear soft coupling, Aq. The fine tuning 
of the electroweak scale is plotted in (a) and (b), the fine tuning of the relic density in 
(e) and (f), and the maximum fine tuning between the two in (c) and (d). The labels 
distinguish between the dominant LSP annihilation mechanisms. The h^ {H^, A^) resonance 
is significant for region 1 (5) with points coloured red (black). Region 3 (4) realises stau 
(stop) co-annihilation with points coloured green (blue), and region 2, coloured purple, has 
increased Higgsino components. 
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Fig 13.91 shows the dependence of fine tunings on the gaugino soft mass. Region 1 
{h^ resonance) fills in a restricted range in mi/2 as already noted from the previous 
plots of neutralino mass. Region 2 (Bino-Higgsino mixing) sets the minimum fine 
tuning for large rn-^^,^, which is slightly raised from the case ignoring dark matter 
constraints. Region 3 (Stau coannihilation) is seen to be restricted to large rn^,2 on 
application of the Higgs mass, which opens a gap in the allowed parameter space 
from 160 to 250 GeV. The discrete lines visible within this region in plots (a) and 
(b) are a result of the discrete steps in tan/3, where larger tan/3 has reduced fine 
tuning. Region 4 (Stop coannihilation) has intermediate Tn-^^,^ ^^^ Region 5 {H^, A^ 
resonance) is on the low fine tuning edge at large m^ /g. 

Fig |3.10l shows the dependence of fine tunings on the trlinear soft parameter. Re- 
gion 1 {h^ resonance) has smallest fine tuning for non-zero Aq, and fills the area just 
above the line of minimum fine tuning when ignoring the dark matter constraints. 
Region 2 (Bino-Higgsino mixing) prefers near zero Aq, with the minimum edges be- 
coming visible on application of the Higgs mass contraint (which doesn't affect this 
region, but excludes points previously overlaying the points of this region). Region 
3 (Stau coannihilation) prefers near zero Aq, but does not extend to the edges found 
without applying dark matter constraints. Region 4 (Stop coannihilation) requires 
large Aq in order to generate the stop mixing as previously argued, where the fine 
tunings are large. Region 5 {H^, A^ resonance) prefers Aq ~ 1 GeV for both fine 
tunings. 

3.5.3 Dependence on LSP composition 

Figs 13.111 to 13.141 plot the fine tunings vs the fraction of the LSP wavefunction given 
by the neutral gaugino and Higgsino states. Fig 13.111 is not helpful for analysing the 
Bino-like regions, but for Region 2 (Bino-Higgsino mixing), the fine tuning of the 
electroweak scale is noticeably reduced as the Bino component is increased. 
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Figure 3.11: Plots of fine tunings vs Bino overlap of LSP. The fine tuning of the electroweak 
scale is plotted in (a) and (b), the fine tuning of the relic density in (e) and (f), and the 
maximum fine tuning between the two in (c) and (d). The labels distinguish between the 
dominant LSP annihilation mechanisms. The /i^ {H^, A^) resonance is significant for region 
1 (5) with points coloured red (black). Region 3 (4) realises stau (stop) co-annihilation 
with points coloured green (blue), and region 2, coloured purple, has increased Higgsino 
components. 
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Figure 3. 12: Plots of fine tunings vs Wino overlap of LSP. The fine tuning of the electroweak 
scale is plotted in (a) and (b), the fine tuning of the relic density in (e) and (f), and the 
maximum fine tuning between the two in (c) and (d). The labels distinguish between the 
dominant LSP annihilation mechanisms. The h? {H^, A^) resonance is significant for region 
1 (5) with points coloured red (black). Region 3 (4) realises stau (stop) co-annihilation 
with points coloured green (blue), and region 2, coloured purple, has increased Higgsino 
components. 
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Figure 3.13: Plots of fine tunings vs Higgsino (neutral down-type) overlap of LSP. The fine 
tuning of the electroweak scale is plotted in (a) and (b), the fine tuning of the relic density 
in (e) and (f), and the maximum fine tuning between the two in (c) and (d). The labels 
distinguish between the dominant LSP annihilation mechanisms. The h^ {H^, A^) resonance 
is significant for region 1 (5) with points coloured red (black). Region 3 (4) realises stau 
(stop) co-annihilation with points coloured green (blue), and region 2, coloured purple, has 
increased Higgsino components. 



109 




\<X°\h2>\' 

(a) Higgs mass unrestricted 




(b) m^> 114.4 GeV 




\<X°\h2>\' 

(c) Higgs mass unrestricted 



1000 




10"^ 10"^ 0.001 0.01 

\<X°\h2>\^ 

(e) Higgs mass unrestricted 




\<x"\h2>\' 

(d) m;^> 114.4 GeV 






100 



A" 




\<x"\h2>r 

(f) mf^> 114.4 GeV 



Figure 3.14: Plots of fine tunings vs Higgsino (neutral up-type) overlap of LSP. The fine 
tuning of the electroweak scale is plotted in (a) and (b) , the fine tuning of the relic density 
in (e) and (f), and the maximum fine tuning between the two in (c) and (d). The labels 
distinguish between the dominant LSP annihilation mechanisms. The h^ {H^, A^) resonance 
is significant for region 1 (5) with points coloured red (black). Region 3 (4) realises stau 
(stop) co-annihilation with points coloured green (blue), and region 2, coloured purple, has 
increased Higgsino components. 
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This follows since for the points of smallest Bino component, the fn-^^,^ is a factor 
of two larger than at the minimum of fine tuning. The heavier spectrum leads to 
increased fine tuning of the electroweak scale. Note that very natural relic densities 
can be achieved in the limit of small Bino component, even with the Higgs mass 
constraint. 

Fig 13.121 shows a strong dependence on the LSP's Wino composition against fine 
tuning. Fine tuning is reduced for a large Wino component. The strips observed in 
Region 3 (Stau coannihilation) is an artefact of the stepped tan/3 scan, with low tan/3 
having larger fine tuning. The fine tuning dependence is similarly observed in Fig l3.13l 
and 13.141 for the Higgsino components. For Region 1 (/i° resonance), increased non- 
Bino components indicates moving into the focus point region, where the fine tuning 
is reached. For Regions 3 (Stau coannihilation), 4 (Stop coannihilation) and 5 {H^, A^ 
resonance), the smaller non-Bino components occur for large ttl^i^, where mixing is 
reduced. The heavy spectrum then leads to increased fine tunings. 

3.6 Summary 

The generation of a dark matter candidate for the mSUGRA model with electroweak 
interactions requires a consideration of the constraints from dark matter physics on 
the allowed parameter space of the model. As the dark matter candidate is mostly 
Bino-like in the majority of parameter space, with small annihilation cross sections, 
the requirement for a thermal relic density in agreement with the WMAP observations 
for the non-baryonic energy density strongly restricts the allowed parameter space. 

The regions of remaining parameter space where increased annihilation cross sec- 
tions were found could be distinguished according to the mechanism present that 
strengthened the annihilation. The mechanisms included resonance via the Higgs 
fields, h^, H^, A^, coannihilation with the staus and stops and residing in the focus 
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point region where significant Bino-Higgsino mixing is present. 

The CDMS II constraint was found to further limit the allowed parameter, being 
most sensitive to the focus point region. The composition of the LSP was inspected, 
from which the phenomenology of the dark matter could be understood. The allowed 
mSUGRA parameter space was presented, demanding consistency with present day 
experimental limits, including the dark matter constraints, for a scan of all mSUGRA 
parameters. 

The fine tunings were plotted and compared for determining the worst fine tuning 
of observables present in the model that are connected to the electroweak physics. 
The CDMS II constraint was found to be only significant below the LEP II Higgs 
limit. The minimum fine tuning with dark matter constraints was 1 part in 15 (for 
which relic density does not saturate the WMAP result and annihilation is via the 
h^ resonance), with rrih = 114.5 GeV. This is only a factor of two larger than that 
found when ignoring the dark matter constraints. In the case of requiring saturation 
of the non-baryonic matter by the neutralino LSP, and the Higgs mass constraint, 
the minimum fine tuning is 1 part in 20 with m^ = 115 GeV (again corresponding to 
annihilating close to the /i° resonance). For points in the focus point, the minimum 
fine tuning of the electroweak scale is 1 part in 30, but for this point the relic density 
fine tuning is 1 part in 40, with m^ = 118 GeV. These results were found without 
application of the Higgs mass bound. 

The dark matter constraints have slightly raised the minimum fine tuning of the 
electroweak scale present, but the naturalness limits (A < 100) are not significantly 
changed for the allowed range for the mSUGRA parameters. The consideration of 
the fine tuning of the relic density is also frequently competitive with that of the 
electroweak scale for the lowest fine tuned points. 

In the near future, different sections of the mSUGRA parameter space will be 
tested by direct dark matter searches and searches at the LHC. In the first run of 
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the LHC which is due to end during winter 2011, an integrated luminosity of lfb~^ 
will be collected at the centre of mass energy of ^/s = 7 TeV. The gluino mass reach 
for the mSUGRA model during this period will extend to 620 GeV |59j. This will be 
able to test the region where the LSP annihilates via the h^ resonance (region 1) for 
which the spin independent cross for scattering of the LSP with a nucleon is up to 
five orders of magnitude below the current experimental limits. The gluino mass for 
region 1 points is 0(400 GeV), but the remaining labelled regions all have too heavy 
gluino masses for detection in the first run. 

For discovery of the gluino, the kinematic edges observed in the distributions of 
its decay products such as charginos and neutralinos and the constrained ratios of 
the mSUGRA spectrum states can be used to demonstrate detection. However, the 
backgrounds are usually significant for the region 1 points and so rapid observation is 
not expected. A useful signal for heavier spectra is the presence of a strong dilepton 
peak from subsequent electroweak boson decay, but within region 1, the mass differ- 
ence of kinematically accessible neutralinos and charginos is insufficient to produce 
on-shell electroweak bosons from cascade decays of the gluino. 

In comparison, the direct dark matter detection experiments probe the region of 
mSUGRA parameter space close to the focus point where there is significant Bino- 
Higgsino mixing in the LSP (region 2). The points in this region have a gluino mass 
in excess of 750 GeV. The points with smallest fine tuning that satisfy the Higgs 
mass bound populate regions 1 and 2, and so these two complementary types of 
experiment provide separate tests of naturalness and either may be expected to soon 
discover natural SUSY physics. 
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Chapter 4 



Improving naturalness through 
non-minimality 



The previous chapters have analysed the degree of fine tuning present in minimal 
supergravity with universal soft mass parameters. This was motivated to produce up 
to date and more precise results over a more complete parameter space scan than the 
analyses previously published. In addition, the minimal scenario was selected for its 
greater transparency in what features control the fine tuning measure. The motiva- 
tions for low energy supersymmetry presented in the introduction did not however 
require a minimality condition. For an appreciation of how naturalness may be mod- 
ified by non-minimal models, this chapter explores such scenarios. In particular, the 
features of models that allow smaller fine tuning with respect to mSUGRA will be 
identified. 

One structure alluded to in previous discussion that may permit reduced fine 
tuning is alternative patterns in the soft mass parameters. For example, consider the 
separate contributions to the electroweak scale from the different soft terms in the 
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MSSM. For tan /3 = 2.5, 



m|/2 = -0.87 l/iol V 3.6 M|-0. 12 M2 + 0.007 M^ 
+ 0.25 M2 M3 + 0.03 Ml M3 + 0.007 M^ Mg 

- 0.71 m^^ + 0.19 ml^ + 0.48 (m^ + m^) 

- 0.34 A^ M3 - 0.07 A^ M^ - 0.01 A^ M^ + 0.09 A^ (4.1) 

where all parameters in the above equation are defined at the UV scale of the theory. 
If all these parameters are independent, then fine tuning will be found to be large. 
However, if say the gaugino masses are fixed according to certain ratios which suppress 
their overall contribution to the electroweak scale, a gaugino focus point is found 
[631 El]. If a non-singlet unified gauge representation is responsible for SUSY breaking, 
non-universal but fixed ratio gaugino masses are generated that may lead to such a 
gaugino focus point. 

An alternative approach for finding improved naturalness is through non-minimal 
particle content. This may relax the experimental constraints, and so open up less 
fine tuned points, help increase the annihilation cross section of the dark matter (or 
provide an alternative candidate) and/or increase the effective quartic coupling of the 
Higgs potential. The investigation of the analytical fine tuning results identified that 
increasing the effective quartic coupling significantly improves the fine tuning. 

To improve naturalness with extra physics, complete model extensions can be 
considered, or if there is a scale separation between the MSSM and the new physics, 
an effective theory can be examined. The latter investigation is now considered for 
extensions of mSUGRA. Once the most efficient couplings for the higher dimension 
operators are ascertained, the models of new physics that can generate such operators 
will be inferred as being natural extensions. 

An understanding of extensions of SUSY models that improve naturalness may 
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become particularly useful if restrictive experimental limits force benchmark models 
such as mSUGRA into increasingly fine tuned regions (for example, if an exclusion 
limit of rrif^ > 121 GeV is found). The premise of natural supersymmetry would then 
suggest specific new degrees of freedom are necessary to preserve the viability of the 
principle of naturalness. 

For model building beyond the MSSM, there are many directions that can be 
followed. The next-to-minimal supersymmetric model (NMSSM) [65] has one extra 
chiral singlet, but an arbitrary number of superfields could be added with appropriate 
quantum numbers and Lagrangian parameters to survive experimental constraints. 
In the following analysis, the model-independent approach of using an effective theory 
with higher dimensional operators [011 EZl ESI ESI EOl ED E2] is adopted in order to 
avoid needing to specify a model. 

4.1 Higher dimensional operators of the MSSM 
Higgs sector 

In this section the effective operators of dimension d = 5,6 are listed that can be 
present in the Higgs sector consistent with the symmetries of the MSSM. These 
operators parametrise new physics beyond the MSSM and affect the Higgs scalar 
potential. Therefore they also affect the amount of fine tuning of the electroweak 
scale, as discussed in detail in the next section. The (i?— parity conserving) d = 5 
operators in the MSSM Higgs sector are: 

A = ^ld'eX{S)iH,H,)^ (4.2) 

C2 = ^f d^e A{S, S^)D° B{S, S^) H2 e-^' D^ C{S, S^) e^' Hi + h.c.(4.3) 
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where S is the spurion field, S = OOttiq, A{S, S'^), B{S, S''), C{S, S'^) are polynomials 
in S, S*^ and mo is the soft scalar mass term in the visible sector. As will be demon- 
strated in Section U2] the first operator can be generated, for example, by integrating 
out massive gauge singlets or SU(2) triplets, while the second is easily generated by 
integrating out a pair of massive Higgs doublets ^71j, all of mass of order M*. 

In [tH [72] it was shown that by using general field redefinitions one can remove 
£.2 from the action. The effect of this is an overall renormalisation of the soft terms 
and of the ii term. Since the fine tuning measure includes the fine tuning with respect 
to each of these soft operators separately, adding £2 cannot reduce the overall fine 
tuning. For this reason only £1 will be included in the discussion of fine tuning with 
(i = 5 operators. 

The £1 term generates the following contributions to the Higgs scalar potential 
on integrating out the auxiliary fields, 

Vi = CiKK{\K\^ + \h2\^) + \C2{Kh2f + ^-C- (4.4) 

where C,^ = 2A/iQ/M„, and C2 = —^X'ttiq/M^ given that X{S) = A + X'S. In the two 
Higgs doublet notation, A5 = Ci and Ag 7 = ^2- At tree level in the MSSM, A5 g 7 = 0. 
There are also d = 6 operators that can be present in addition to the MSSM 
Higgs sector. These are suppressed relative to the d = 5 operators by the factor 
1/M^,. However they may give contributions to the Higgs potential enhanced by 
tan/3 relative to the d = 5 so cannot be ignored at very large tan/3. The list oi d = 6 
operators is (see also [HZIEH]): 

Oi = -^J d^'e Z,{S,S^){Hle^^Hi)\ ^ = 1,2. 

O3 = ^jd'eZsiS,S^)iHle''^H,)iHle'''H2), (4.5) 
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The above operators can be generated by integrating a massive f/(l) gauge boson or 
a SU{2) triplet. 



04 = 


1 


o, = 


1 


Oe - 


1 
Ml 


O7 = 


1 


Os = 


1 
M2 



d^e Z^{S, S^) {Hi e^i Hi) {H2 Hi + h.c.) 
d^e Zq{S, S^) {Hi e^^ H2) {H2 Hi + h.c.) 
d'^e Z^{S, 0) W" Wo, {H2 Hi) + h.c, 



d^e 



Z^{Q,S^) {H2Hi)' + h.c. 



(4.6) 



where Wa is the supersymmetric field strength of a vector superfield of the SM gauge 
group. As an example, O4 can be generated by integrating a gauge singlet. 



O9 

Oio 
On 

012 



1 

MS 



d^e Zg{S,S^) Hlv\^'V^Hi 



-^ J d'e Zio(5, S^) Hi V' e^^ V^ H2 
-^ j d^e Zii{S, S^) Hi e^i V^ 1^"^ ffi 
/" d^e Zi2{S, S^) Hi e^^ V^ iy« i/2 



1 

Af2 



(4.7) 



where V° acts on everything to the right and V^ H^ = e~ ' Da e^H^. In addition to 
the spurion dependence in the wavefunctions Zi{S, S^), extra {S, S^) dependence (not 
shown) can be present under each derivative Vo in eq. (l4.7p . in order to ensure the most 
general supersymmetry breaking contribution associated to these operators. One may 
use the equations of motion to replace the operators involving extra derivatives by 
non-derivative one^j. Note that when computing the fine tuning measure eliminating 



^ Setting higher derivative operators onsheU is a subtle issue in this case. One can also use general 
spurion-dependent field redefinitions to "gauge away" (some of) these operators, using the method 

of imizs]- 



118 



a particular operator will lead to correlations between the remaining operators that, 
strictly, should be taken into account. 

As the determination of fine tuning depends on using RGEs between the gauge 
unification scale and the electroweak scale, one cannot consistently consider all di- 
mension six operators when M^ is between these scales. However, this is the case 
of interest in order to find significant effects on electroweak physics. This limitation 
arises because the physical degrees of freedom that have been integrated out to gener- 
ate these operators must be accounted for in the running when above their threshold 
scales. For this reason, numerical results of fine tuning in cases where dimension 
six operators are important are not presented. As states that do not significantly 
affect the running of MSSM states, such as gauge singlets, can generate dimension 
five operators, reliable fine tuning results in the presence of dimension five operators 
can still be determined. 

4.1.1 Validity of effective theory 

The operator analysis used here has a limited range of validity because it corresponds 
to integrating out new heavy degrees of freedom. If the mass of these degrees of 
freedom is not much above the energies being probed, the operator analysis breaks 
down and one must deal with the new degrees of freedom directly. The mass, M^,, at 
which this happens corresponds to the point where high dimension operators are not 
suppressed relative to low dimension operators. A measure of this may be obtained 
by dimensional analysis in which the operator matrix elements are taken to be deter- 
mined by the energy scale being probed. Applied here, this implies that the operator 
analysis is reliable provided mo, // ^ M^. 

A potential fault in this estimate of the range of convergence occurs because higher 
dimension operators may have anomalously large matrix elements. An example of 
this occurs for the dimension-six operators listed in Section 14.11 Consider the first 
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dimension- six operator in eq.f l4.5p 



<^^^lk 



S^S{Hle''H,y\^D^-]^\h,\' (4. 



for Ci ~ 0(1). This should be compared to the leading quartic Higgs term coming 
from the dimension-five operators in eq (14 ■2p that contributes at O ijp- |/i2| /ii/i2)- 
One may see that the relative magnitude of the dimension-six to dimension-five contri- 
butions is O i 2~M~ ^^^ (^ ) ■ Thus, strictly, the region of validity of the dimension-five 

2 

operator analysis is ^ "^^/^ tan /3 ^ 1. However, as discussed in the next section, the 
new physics generating this dimension-six operator is different from that generating 
the dimension-five operator and so their coefficients should be uncorrelated. In this 
case the addition of higher dimension operators can reduce the fine tuning for some 
region in parameter space so that the analysis with dimension-five operators only 
will provide a useful upper bound even in regions where dimension-six contributions 
are significant. For this reason the region of validity of the dimension-five operators 
analysis is better described by the original mo/M* <^ 1 and fio/M^. <^ 1 condition. 
This keeps the corrections coming from operators with correlated coefficients small. 
In the following quantitative analysis, only the effects of dimension-five operators 
are included. The convergence criterion found above gave mo/M^: ^ 1 and ^o/M^, ^ 
1. These bounds are comfortably satisfied in the following when the restriction, 
rrio/M^, iiq/M^ < 0.035 is made, giving upper values Ci,2 < 0.07. 

4.1.2 Fine tuning results from parameter space scan 

The fine tuning of the electroweak scale is considered here, using the analytic formulae 
of Chapter 2 for the fine tuning measure, Higgs potential and associated spectrum. 
The dark matter physics in the presence of higher dimension operators is not consid- 
ered as the increased set of relevant dimension five operators for LSP annihilation and 
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Figure 4.1: Left figure (a): the MSSM fine tuning A as a function of rrih', Riglit 
figure (b): the fine tuning in the MSSM with d = 5 operators in terms of rrih, with 
Ci = C2 = 0.03. In both figures, the top pole mass considered is rrit = 174 GeV for blue 
(dark blue) areas and rrit = 171.2 for yellow (red) areas, respectively. Larger m^ input 
(blue) shifts the plots towards higher rrih by 2-5 GeV. In both figures the parameters 
space scanned is: 1.5 < tan/3 < 10, 50 GeV < mo,mi2 <1 TeV, —10 < Aq/ttIq < 10, 
without application of the experimental constraints. 

LSP-nucleon scattering significantly modifies the allowed parameter space for which 
an analysis is beyond the scope of this discussion. One of the key results of this sec- 
tion is a demonstration of how the dependence of fine tuning tuning on Higgs mass 
is modified by non-minimal physics. In the mSUGRA case, Higgs masses of greater 
than 121 GeV required A > 100. This value is quite close to the current Higgs mass 
bound. The following results indicate that non-minimal physics permit larger Higgs 
masses with smaller fine tunings. 

In a similar approach to the previous results, a parameter space scan is performed 
and the resulting fine tuning values plotted for the points tested. Due to the interest in 
identifying relative changes and general features introduced by higher dimension op- 
erators, the experimental constraints are ignored here. This also avoids complications 
that extra higher dimension operators would affect the experimental constraints. 

To ensure that the calculations are appropriate (ie neglected contributions are 
insignificant at the level of precision considered), the range of tan/3 is restricted to 
low values between 1.5 and 10. The Higgs potential is also only evaluated at one 
loop leading log order. The parameter space scan is further defined by, 50 GeV < 
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Figure 4.2: Left figure (a): the fine tuning A as a function of ruh- A of MSSM is 
plotted in light blue (mj = 174 GeV) with an orange edge (mj = 171.2); A of MSSM 
with d=5 operators with Ci = C2 = 0.03 is plotted in dark blue {rrit = 174 GeV) with a 
red edge (m^ = 171.2). Right figure (b): similar to figure (a) but with Ci = C2 = 0.05. 
Non-zero or larger Q (dark blue and red areas) shift the plots to higher rrih to allow 
a reduced A for higher rrih- In both figures 1.5 < tan/3 < 10; 50 GeV < mo, mu < 1 
TeV, —10 < Aq/itiq < 10, without application of the experimental constraints. 

7710,17112 < 1 TeV, and the trilinear coupling is varied according to, —10 < Aq/tuq < 
10. The signs of Ci 2 ^^^ always chosen such that fine tuning is reduced with respect 
to the mSUGRA case. 

The results are shown in Figures 14.11 to 14.31 Note that in these figures the struc- 
ture apparent at small A and large rrih is a scanning artefact due to discrete steps in 
parameter space. As a benchmark. Figure I4.ir a) shows the fine tuning of the elec- 
troweak scale. A, in mSUGRA given the method used in this section. It was observed 
in Chapter 2 that the right hand edge of this plot is insensitive to the experimental 
constraints. The same dependence of fine tuning on Higgs mass is then reproduced. 

Figure 14.1( b) shows A for the case of the MSSM with dimension-five operators 
added, with Ci = C2 = 0.03. The dominant effects in Figure I4.ir b) are mostly due 
to the effect of non-zero (i, which comes from the supersymmetric part of the higher 
dimensional operator. One may see a systematic shift of the allowed region to higher 
TJih which (for positive Q) is driven by an increase in the quartic Higgs coupling which 
appears in the denominator of the fine tuning measure. The overall result is that the 
minimum amount of fine-tuning A in the presence oi d = 5 effective operators is small. 
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Figure 4.3: As in Fig J4.2l but with: left figure (a): Ci = 0.07, ^2 = 0; right figure (b): 

Ci = 0, C2 = 0.1 

of order A ^ 6, for rrih from 95 to 119 GeV. Therefore non-zero Q can accommodate 
larger rrih while keeping a A significantly smaller than in the MSSM, for a similar 
value of rrih. To illustrate the change more directly both plots are superposed in 
Figure H2]^a). The effect is enhanced for larger operator coefficients, as may be seen 
in Figure U2Ib), where A is presented for Ci = C2 = 0.05, again shown relative to the 
MSSM case. One can see that in this case, values of rrih > 114.4 GeV can have a low 
A ^ 6. Therefore A can be significantly reduced from the MSSM case, for a similar 
THh- This conclusion is further supported by the plots in Figure H^ where other values 
for Q are considered. From all plots shown, one sees that A < 10 is easily satisfied 
for values of the Higgs mass that can be as large as 130 GeV, depending on the exact 
values of Ci- 

Note that, in the MSSM, A increases for low tan/3 (^ 10) and rrih above the 
LEPII bound. However, the effect of the d = 5 operators is important for low tan /3 
and in their presence A actually decreases for low tan/3. Thus the reduction in the 
fine tuning at very low tan/3 relative to the MSSM case is much more marked than 
that shown. 

The lower amount of fine tuning in the presence of effective operators is due to 
two effects. The first, already mentioned, is the presence of additional quartic Higgs 
couplings enhancing the denominator which determines the Higgs via f ^ = —m?/\ 
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thus allowing for a smaller electroweak breaking scale. The second is the fact that 
higher dimensional operators add a tree level contribution to the Higgs mass, which 
reduces the need for large quantum contributions, and therefore the fine tuning. 

The scale of new physics that is generating this increase in Higgs mass and reduc- 
tion of fine tuning is given by, 

M* ^ 2/io/Ci ~ (40 to 65) X /io, 0,2 = 0.05 to 0.03 (4.9) 

For yUo between the electroweak scale and 1 TeV, this shows that large values of M* 
are allowed: M, ^ (5.2 to 8.45) TeV for /io = 130 GeV and M, ^ (8 to 13) TeV for 
Hq = 200 GeV. For larger /io one obtains values of M^, above the LHC reach. Finally, 
for Ci = 0.07 but with (2 = 0, one has M=, ^ 30 x /io and A < 10 for nih ^ 130 
GeV. Thus, the EW fine tuning is small A < 10 for 114 < nih < 130 GeV, for rather 
conservative values of 0,2- To relax these values, a rough rule is that an increase of 
(i by 0.01 increases nih by 2 to 4 GeV for the same A. 

4.2 The origin of "new physics" 

The presence of a higher dimension operator signals new physics and it is important to 
ask what this new physics can be. In the context of new renormalisable interactions 
it may come from the effects of new chiral superfields or from new gauge vector 
superfields. Consider chiral superfields first. One may readily obtain the d = 5 
operator of eg. (14. 2 p by integrating out a gauge singlet or a triplet [70]. Consider the 
case of a massive gauge singlet X with Lagrangian 

Cx= fd^eX^X+l f d^e fiHiH2 + K X H1H2 + I M, X^ +h.c.\. 
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For M^, 3> /i, rriQ, one may use the equation of motion to integrate out X, giving, 



D^X^ + MX + W' + XW" = 



(4.10) 



where W = X^XH^H^ + fiH^H^ and W = dW/dX, and so on for W" (which is zero 
in the above example). A series expansion of this equation of motion generates, 



X 



w 



M 4\M\ 



^^+o(i9.) 



(4.11) 



which can be used to substitute into the original Lagrangian. This generates dimen- 
sion five and six operators up to 0{M^^) terms. 



^•effective 



M2 



d^e\H^H^f + 



^ J d'e {H,H,f + h.c. 



(4.12) 



The supersymmetry breaking terms associated with this operator are obtained by 
replacing A -^ X{S) giving the d = 5 operator of interest. Note that Cx has a similar 
form to that of the NMSSM. However, in the NMSSM, the singlet field has mass of 
order the electroweak breaking scale and cannot be integrated out whereas here, the 
singlet mass is assumed to be much larger than the EW scale. 

However, the origin of the d = 5 operator cannot be uniquely ascribed to a gauge 
singlet field. Indeed it may equally well point to the existence of SU{2) triplets 
[73| [HI [75], Ti 2,3, of hypercharge ±1, 0. In this case a Lagrangian of the form 



C, 



dW 



-i^v 



UVr 



TyTi+Tle''T2 



d'e 



fiHiH2+M,TiT2+XiHiTiHi+\2H2T2H2 



+h.c 



which generates an equivalent dimension five operator as in eq (I4.12p except that A^ 
is replaced by A1A2. More generally, one can generate the d = 5 operator through a 
combination of both gauge singlets and triplets. However note that the pure singlet X 
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case has the advantage of not affecting the gauge couphngs unification (at one- loop), 
which is not true for the SU{2) triplet. 

What about additional, massive, SU{2) doublets that couple to the MSSM Higgs 
sector? One may readily show that integrating them out does not generate, to lowest 
order in 1/M*, a dimension five operator of the type in eq f l4.f2p . 

There remains the possibility that the new physics is due to the effect of new 
massive vector gauge superfields. The simplest example is the case where there is 
a new ^(1)' gauge symmetry under which the Higgs sector is charged. This brings 
extra quartic contributions to the scalar potential that are expected to reduce the 
fine-tuning [761 [771 [78] . Assuming the f/(l)' is broken at M^: one obtains the effective 
Lagrangian to leading order in inverse powers of M^, given by 



/-effective _ _gj_ / Mn 
'~'U{1) ~ A/2 ' " ^ 



'2 I - r ... ... i2 



q,Hle''H, + q,Hle''H, 



where g' is the U{iy coupling and gi_2 are the charges of the Higgses under f/(l)' 
(?! + ^2 = 0). Note that, after including the associated supersymmetry breaking 
operators, this corresponds to the d = 6 effective operators [70j of eq. (14. 5 p and that 
no d = 5 operators are generated. 

In summary, the requirement that the SUSY extension of the MSSM should not 
have significant fine tuning may indicate the presence of the d = 5 operator of eq. (14. 2 p 
which, in turn, suggests the presence of a massive gauge singlet and/or a SU{2) triplet. 
This is the simplest interpretation based on new renormalisable interactions but other, 
more complicated possibilities to generate the d = 5 operator may be possible. 
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4.3 Further remarks on fine tuning 

Effective field theory approaches to the fine tuning of the electroweak scale were used 
before in models of low SUSY breaking scale scenarios ^21j where both d = 5 and 
d = 6 operators were included. The model in [2T] introduces supersymmetry breaking 
through coupling of MSSM states to a SM singlet field responsible for supersymmetry 
breaking. After integrating this field out, in addition to the d = 5 operator considered 
here, there are correlated contributions from the d = 6 operators. Using this, the 
authors find the fine tuning can be very small even for an arbitrarily high Higgs 
mass, provided the scale of supersymmetry breaking is less than 500 GeV. 

How does this analysis relate to the one presented here? The examples given 
in [21] are found varying the ratio fh/M in the range 0.05 to 0.8 where fh is the 
supersymmetry breaking scale and M is the messenger mass. For ffi/M small, the 
fine tuning is close to that in the MSSM but reduces rapidly for rh/M large; in 
this latter case the fine tuning actually reduces as the Higgs mass increases. This 
range of values for fh/M corresponds to a choice of our ttiq/M^ and /iq/M^, in a 
similar range. The upper value strongly violates the criterion for applicability of the 
operator analysis argued here and is a factor of ~ 10 larger than the value chosen 
in Figure H^ a). Ignoring, for the moment, the fact that the contributions of higher 
dimension operators are expected to be large for this choice of mediator mass, it can 
be asked what this choice of mediator mass in this analysis would give for the Higgs 
mass consistent with small A. Since the change in the upper bound on the Higgs mass 
roughly scales with the coefficient of the rf = 5 operator, this would allow a Higgs 
mass in the region of 276 GeV, much larger than the earlier conservative estimates 
presented here. However, as has been stressed, for this value of the messenger mass the 
operator analysis breaks down and one should do the analysis including the messenger 
fields explicitly. 
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4.4 Summary 

The LEPII lower bound on the Higgs mass places MSSM Higgs physics at the forefront 
of supersymmetry phenomenology. While this bound can be satisfied by including 
the MSSM quantum corrections, it (re) introduces some amount of fine tuning in the 
model. To reduce the fine tuning may require new physics beyond the MSSM which 
can be parametrised by higher dimensional operators. In this chapter, an effective 
field theory framework with d = 5,6 operators is considered in the MSSM Higgs 
sector, giving a model independent approach to the fine tuning problem. 

Fine tuning proves to be very sensitive to the addition of higher dimensional 
operators and this is mostly due to extra corrections to the quartic couplings of the 
Higgs field. For the case of dimension-five operators it was showed that one can 
maintain a reduced fine-tuning A < 10 for a Higgs mass above the LEPII bound 
and as large as rrih ~ 130 GeV, for the parameter space considered, with low tan/3 
(tan/3 < 10). The scale of new physics M^, responsible for the reduction in fine tuning 
can be rather large, for example M^, ^ 2/io/Ci ~ (40 to 65) x/iq, for ^12 = 0.05 to 0.03, 
and M* ^ 30 x /iq for (i = 0.07, (2 = 0. For values of /iq between the electroweak 
scale and 1 TeV, these results show that large values of M^, are allowed; in the 
former case M^ ^ (5.2 to 8.45) TeV for /xq = 130 GeV and M^ ^ (8 to 13) TeV for 
fj,Q = 200 GeV. For larger /io, larger values of M^, are possible, even above the LHC 
reach. These results follow from rather conservative choices for the coefficients of the 
quartic couplings induced by the dimension-five operators, to ensure the convergence 
of the effective operator expansion. 

The numerical analysis included the effect of dimension-five operators only. These 
give the leading corrections at low tan/3, being proportional to 1/M*. However, 
dimension-six operators, suppressed by itIq/M^ or fi^/M^, give contributions that 
can be enhanced by large tan/3; for (tan /3 mo) > M* or (tan/3/io) > M* these will 
be the leading terms. However, due to what would be an inconsistent application of 
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renormalisation group equations, general dimension six operators cannot be reliably 
tested for fine tuning with this methodology. Depending on the physics that generates 
the dimension five operators, the results determined here may also not be appropriate. 

Of course the crucial question is what is the origin of the physics beyond the 
MSSM giving rise to these operators? The dimension-five operator can be generated 
by a gauge singlet superfield or a SU{2) triplet superfield of mass of 0{M^) coupling 
to the Higgs sector. The dimension-six operators can be generated, for example, by an 
extra gauge symmetry with a massive gauge supermultiplet or additional (Higgs-like) 
SU{2) doublet supermultiplets of mass (9(M^,). 

Returning to alternative methods for reducing naturalness, in addition to non- 
minimal particle content affecting the Higgs potential, experimental constraints may 
be evaded. Alternatively, the pattern of high energy soft SUSY breaking parameters 
may generate special features such as focus points at the electroweak scale. These are 
the main mechanisms for improving naturalness. If the natural parameter space of 
a minimal model is ruled out, such mechanisms should then seriously be considered, 
having initally motivated low energy supersymmetry using a naturalness principle. 



129 



Chapter 5 



Conclusions 



Supersymmetry has been argued to be the most general spacetime symmetry consis- 
tent with a selection of assumed physical principles. It is then appealing to embed a 
model of nature with such a symmetry. Furthermore, there are theoretical motivations 
that lead to expecting the physics of supersymmetry to be realised around the TeV 
scale which is at our current experimental frontier. The motivations include providing 
a natural solution to the hierarchy problem of the electroweak and Planck scales and 
a dark matter candidate which may saturate the inferred non-baryonic matter energy 
density. A thorough understanding of these models may lead to imminent discoveries. 

Having motivated supersymmetry by invoking a principle of naturalness, it is then 
vital to test the allowed regions of parameter space for the presence of a satisfactory 
degree of fine tuning. The parameter space can also have limits placed upon it for 
remaining natural. The phenomenology of the natural regions of parameter space 
can be identified, and search strategies optimised for testing the existence of natural 
theories. 

The work presented here has identified the natural regions of minimal supergravity, 
provided naturalness limits on the parameter space and spectrum and discussed the 
phenomenology associated with the natural realisations. The analysis was completed 
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using the latest experimental constraints, and tested the fine tuning of the electroweak 
scale at the 2 loop order for a scan over all parameters of the theory. The fine tuning 
of the thermal relic density, controlled by the electroweak interactions, were similarly 
tested and found to compete with the former fine tuning in the regions of lowest fine 
tuning. 

A fine tuning of 1 part in 15 was found to be consistent with all current experi- 
mental constraints for the mSUGRA model. This is much smaller than the commonly 
repeated value of 1 part in 100 for MSSM theories. The low fine tuning occurs in 
the region around the scalar focus point, which is a special property obtained for a 
family of structures of the UV soft scalar masses, for which the universality condition 
within mSUGRA belongs. The MSSM spectrum around this low fine tuned region is 
similar to the benchmark point SPS2, for which the phenomenology has been widely 
investigated. 

Methods for improving the naturalness found in mSUGRA through non-minimal 
models was also discussed. These included special patterns in the fundamental input 
parameters that lead to an insensitivity of the electroweak scale to these parameters, 
or non-minimal physics that allows evasion of the experimental constraints. 

Non-minimal physics may also increase the effective quartic coupling of the Higgs 
potential simultaneously increasing the Higgs mass and reducing fine tuning. A model 
independent analysis was presented using an effective theory formalism to demon- 
strate that new physics beyond mSUGRA with an energy scale beyond the design 
reach of the LHC can still allow a Higgs mass of 130 GeV with fine tuning less than 
1 part in 10. 

To conclude, a large group of natural supersymmetric theories remain experimen- 
tally viable which includes minimal supergravity. If the natural parameter space of 
a given model is fully excluded, this will motivate introducing further new physics 
or different UV completions. The phenomenology of natural minimal supergravity 
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suggests that experiments will soon be in a position to discover its physics, allowing 
a new Standard Model to be written. 
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Appendix A 



Renormalisation group equations 



Using the notation from eqs fll.49p and (11.671) . the 1-loop renormahsation group equa- 
tions (RGEs) are given by 



dt 



dt 



9a 

167r 



9a M a 

87r2 



^yijfe 


yijp 


dt 


167r2 


d^/i 
dt 


167r2 


dA^i^ 
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dt 
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dB'J 
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1 
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(A.7) 



133 



where the MSSM state masses are assumed to be neghgible, and a SU{5) gauge 
normahsation is chosen. The SU{5) normahsed hypercharge, Y', is related to the 
hypercharge in Table [TT^ Y, by Y' = YJ^. To keep the Lagrangian invariant, the 
hypercharge gauge coupling must then also be rescaled, gyi = QyX \- 

The group theoretic factors are given by: and for SU{n), the S^ (1) , C„ (1) = 0. 



SU{n) 


U{l)y 


C{G) = 


n 


C{G) = 


C„(n or fi) = 


n2-l 
2n 


C,{Y) = Y' 


S'„(n or n) = 


1 

2n 


SAY) = Y' 



Table A.l: Group theoretic factors for SM symmetries 
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